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Abstract 



It is known that there is a phenomenon called "parameter degeneracy" in neutrino oscillation 
measurement of lepton mixing parameters; A set of the oscillation probabilities, e.g., P(yu — > v e ) 



the values of #13 and 5. With use of the approximate form of the oscillation probability a la 
Cervera et at, a complete analysis of the eightfold parameter degeneracy is presented. We propose 
a unified view of the various types of the degeneracy as invariance of the oscillation probabilities 
under discrete mappings of the mixing parameters. Explicit form of the mapping is obtained either 
by symmetry argument, or by deriving exact analytic expressions of all the degeneracy solutions 
for a given true solution. Due to the one-to-one mapping structure the degeneracy solutions are 
shown to form a network. We extend our analysis into the parameter degeneracy in T- and 
CPT-conjugate measurement as well as to the setup with the golden and the silver channels, 
P{y e — > u^) and P(y e — > v T ). Some characteristic features of the degeneracy solutions in CP- 
conjugate measurement, in particular their energy dependences, are illuminated by utilizing the 
explicit analytic solutions. 
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I. INTRODUCTION 



After establishing the neutrino masses and the lepton flavor mixing pQ by the atmo- 
spheric j2], the solar [3], and the reactor experiments [I], which is further supported by 
the accelerator experiments E], there seems to exist a consensus that the next step is 
to measure #13 and 5, the remaining unknowns in the MNS matrix, and to determine the 
neutrino mass hierarchy. It was proposed that if #13 is relatively large an intense neutrino 
beam from nuclear reactors can be used to measure it by using the near-far two-detector 
setting |8] . Alternatively, or complimentarily, the accelerator search for nonzero #13 has 
advantage of potential possibility of extending it to CP violation search. The reactor [91412"] 
and the accelerator experiments [131 E] are either ongoing or in construction to look for 
effects of nonzero #13. 

It is well known, however, that detection of CP violation due to the lepton Kobayashi- 
Maskawa phase 5 [U], being the genuine three flavor effect, is suppressed by the two small 
factors, the ratio Am^/Am^ [51 161 [16HT8] and the value of #13 bounded from above [T9H22] . 
Therefore, high precision experiments are inevitably required to measure CP violating phase 
5. Once precision measurement becomes the necessity, the experiment is better characterized 
as a simultaneous determination of #13 and 5. It is because even though one enjoys prior crude 
knowledges of magnitude of 6*13 (assuming it relatively large), which certainly propels CP 
measurement, the required precision for detecting tiny effects of 5 necessitates simultaneous 
measurement of #13 in a precision far beyond the previously achieved ones. 

It is known that a set of measurement of the oscillation probabilities, e.g., Piv^ — > v^) and 
its CP conjugate P{v^ — > u e ) at a particular neutrino energy, no matter how accurate, does 
not determine uniquely the values of #13 and S, the problem of parameter degeneracy [2BT125] . 
The nature of the degeneracy can be understood as the so called intrinsic degeneracy [23] 
duplicated by the unknown sign of Am| x [21] and 823 octant [25], which entails the total 
eightfold degeneracy if 6*23 7^ tt/4. The feature can be seen in Fig. [TJ Some features of the 
degeneracy were further discussed in [2TH29] . whose first two noticed special features that 
appear in the vacuum oscillation maximum. The notorious feature of the degeneracy is that 
difference between the true and the fake solutions can be so small that their distinction is 
extremely difficult, rendering resolution of completely different physical pictures, e.g., the 
mass hierarchies, untenable. Or, in the other cases, the difference between the true and the 
fake values of 5 is so large to confuse CP violation with CP conservation. 

It is the purpose of this paper to give a complete analysis of the parameter degeneracy 
in neutrino oscillations. To achieve a unified understanding of the phenomenon, we present 
and advance a new view of the degeneracy as an invariance of the oscillation probabilities 
under discrete mappings of the mixing parameters. With use of the approximate form 
of the oscillation probability obtained in [30] we present a self-contained derivation (re- 
derivation in CP conjugate case) of the analytic expressions of all the degeneracy solutions 
as functions of the true solutions, which supplies the explicit form of the mapping. Having 
the analytic solutions of the eightfold degeneracy with the proper convention at hand, we 
demonstrate that they form a solution network, the one-to-one correspondence structure 
between solutions, which will be pictorially represented in Fig. [4] in Sec. [V] 

We use the analytic solutions to make plots of the differences between the true and the 
fake solutions to illuminate the global overview of the degeneracy. We illuminate, by using 
the plots, the characteristic features of the degeneracy and reveal the reasons why and how 
the sign-Am^ and the 623 octant degeneracies are robust against the spectrum analysis. 
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FIG. 1: Left panel: An illustrative example of the eightfold degeneracy is represented pictorially 
(as first appeared in [26J) in terms of the bi-probability plot in P — P CP space [21] ■ Right panel: 
The parameters (sin 2 2^13, 5) of the true solutions and the clone ones II— VIII are presented as 
numbers and also by the dots in sin 2 29n — 5 space. The correspondence between the ellipses and 
the solution labels are made manifest by using the same color lines/symbols in both panels. 



We note that the first attempt toward analytic solutions of the degeneracy was pursuit 
by the authors of [27] who obtained the solution for the intrinsic degeneracy. Then, the 
similar analysis was extended in [3T] to include the degeneracy solutions which involve the 
Amjj-sign and/or the 6*23 octant flips. 

One may ask; Why is the parameter degeneracy defined as above way so relevant? Math- 
ematically speaking, the degeneracy is easy to solve; Repeating measurement at one more 
energy (or baseline), or adding a different oscillation channel immediately solves degeneracy. 
Or, if the spectrum information is available it may be more powerful to resolve the degen- 
eracy. Despite these valid reasonings, unfortunately, at least some type of the degeneracy is 
shown to be robust and survives in varying experimental settings. One of the reasons for it 
is, as we will see in Sec. VI, that the energy dependence of difference between the true and 
the fake solutions is so mild that spectrum information is not powerful enough to resolve 
the degeneracy. 

We emphasize that need for resolution of the degeneracy is not only because precision 
measurement is always desirable, but also because, far more importantly, it leaves e.g., the 
neutrino mass hierarchy undetermined after huge experimental efforts. We want to warn the 
readers that we will not try to discuss how the degeneracy can be lifted by assuming concrete 
experimental settings. 1 Rather, we focus in on a complete understanding of structure of the 
degeneracy. The explicit analytic solutions and knowledges of structure of the degeneracy 
should serve for clearer understanding of the experimental data taken in precision measure- 
ment in the future. It will be definitely called for if future neutrino experiments reveal 
features that may not fit in into the standard three-neutrino mixing to clearly discriminate 



The analysis presented in this paper may be regarded as "pathological analysis of neutrino oscillation" . 
It would not tell us directly the experimental method for solving the degeneracy, but as in the case of 
pathology of human body, the clearer understanding of the disease may ultimately provide with us the 
correct recipe for its resoluion. 
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confusion by the degeneracy from new effects beyond the standard three-flavor mixing. For 
example, if they are so powerful to detect neutrino's nonstandard interactions (NSI) [32H36] 
(see |37J for further references therein), the event structure will be modified by the new 
ingredients and enriched with new type of the degeneracies [35J E2] • For the importance of 
lifting the degeneracy, quite naturally, a great amount of efforts were devoted to investigate 
how it can be done. The references [4*01464"] are nothing but only a small subset of them. 

In the next section, we define our machinery and introduce a new view of the degeneracy 
as invariance under the discrete transformations. Then, in the following two sections (III 



and IV), we first give a complete treatment of the parameter degeneracy with CP-conjugate 
measurement. In Sec. |Vj we present the explicit form of the discrete mapping and com- 
plete our understanding of the structure of the degeneracy. It will be supplemented by the 
discussions of degeneracy in T-conjugate (Sec. VII), the golden-silver (Sec. VIII), and CPT- 
conjugate (Sec. IX) measurement combining the v e — >■ (golden) and the u e — > v T (silver) 
channels. 



II. NEUTRINO OSCILLATION PROBABILITY AND ITS INVARIANCE 

In our analysis in this paper, we rely on the approximate formula for the appearance 
oscillation probability derived by Cervera et al. [50] ■ A simple way of deriving the formula 

is to use perturbative framework in which si 3 is assumed to be of order = e and 

keep the terms up to second order in e. For a review of this method, see e.g., [38]. Here, 
Amjj = m| — mf = 1,2,3). However, in this paper, we take an attitude to utilize 
the formula as far as it is reasonably accurate, even outside the region of validity of the 
perturbative ansatz. In fact, it is known that the formula gives a reasonable description of 
the oscillation probability even with larger values of s i3 [65J. 



A. Approximate formula of the neutrino oscillation probability 

To present the formula for the v e appearance probability in a compact way we use the 
simplified notations. We summarize them together with their magnitudes for convenience 
of the readers: 



s = s 13 



A;* = 



AmjiL 



\Amj 



2 



(t) ( JLY 1 



AE 

A = ^ — 0.27 ( — 3^ ( — — — ^ , (1) 

4£ V2.8g/cmV V 1000km ^ 1 ; 

where a = 2y/2GpN e E, the well known coefficient related to the index of refraction of 
neutrinos in matter [32]. Gf is the Fermi constant, p and N e = p/m^ with being the 
nucleon mass denote, respectively, the averaged matter and the electron number densities 
along the neutrino trajectory in the earth, and we have assumed that the electron fraction 
Y e is 0.5. By using the definition of as positive definite quantities we choose to display 
explicitly the sign of Arwgj as ± signs (sometimes called as the hierarchy signs) in the 
equations. 
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It may be useful to remember the ratio between the vacuum and the matter parameters 



for understanding the feature of the degeneracy solutions in Sec. [VT 
A 



0.085 



^31 



|Am 2 J 



2.5 x 10" 3 eV 2 



2.8g/cmV \lGeVj ' 



(2) 



Therefore, in typical low-energy superbeam [661 - 168"] experiments the ratio is small, A/A31 ^ 
0.1, whereas in neutrino factory [691 ITU] with baseline of several thousand kilometers the 
ratio is large, A/ A31 ~ 3 — 6. 

The oscillation probabilities of the neutrino flavor conversion processes — > v e , its CP- 
conjugate channel — » u ei the T-conjugate channel v e — > z/^, and the CPT-conjugate 
channel v e — > z/ M , in matter are given under the constant matter density approximation as 
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(6) 



where ± indicates the mass hierarchy, namely, the normal and the inverted mass hierarchies 
for the positive and the negative Am^, respectively. The functions X±, Y±, and Z are 
defined by 



X± 
Y ± 
Z 



4s 



2:-! 



A 3 i sin(A 31 =f A) 



1 2 



(A31 =F A) 
±2^/X ± P e = ±4sin26 l 1 2C23S23 



A 3 i sin(A 31 t A) 
(A 3 i =F A) 



A21 sin A 
A 



^3 sin 2 26*12 



A21 sin^ 
A~ 



1 2 



(7) 



Their forms imply that the oscillation probability can be written as P = \y/Xs + 
e i(<5±A 3 i)y^|2^ ^jch allows simple interpretation of the (^-sensitive term as an interference 
between the atmospheric and the solar scale oscillations. X and Y in Q and ^ are related 
to X and Y as 



Xfc(a) 
Y±(a) 



X±( 
Y ± ( 



-a) 
-a) 



X T {a) 
-Y T (a). 



(8) 



In our discussions in this paper, it is crucial to note the relation [29] between the coefficients 
X± and Y±: 

Y+ Y_ 

(9) 



which follows from the definitions. Notice that (J8|) means that the same relation as ^ holds 
also for X and Y. 
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In this paper, our emphasis is placed on the oscillation channels between and v e 
and their anti-particles. To have a clearer view of the structure of parameter degeneracy, 
however, we will include the u e — > v T appearance channel, which is sometimes called the 
"silver channel" [71]. See Sec. VIII The oscillation probability P{v e — > v T ) is given by 



P b = P(v P . ->• v T 



cot 2 2?> X±s 2 - Y±s cos (5 =F A 3i ) + tan 2 9 23 Z. 



(10) 



B. Parameter degeneracy as an invariance of the oscillation probability under 
discrete mapping 



It is not so well recognized that the appearance oscillation probability in matter under 
the Cervera et al. approximation has an invariance. Namely, it is easy to show that the 
oscillation probabilities P, P T , and P s defined in rf3|) , (pi) and ( 10 ), respectively, with positive 



Am 2 ! is invariant under the transformation 



-Am 2 1; 




(11) 

Notice that the first transformation transforms X + and Y + into X_ and Y-, respectively, 



and use has been made of the key relation J9 



Similarly, the oscillation probability with 



negative Am| x has an invariance similar to (11), replacing the second one with s — > , ,, v 

It is nothing but a generalization of the invariance of the oscillation probability in vacuum to 
that in matter, which was used to show the existence of degeneracy solutions with differing 
sign of Am 



31 



Then, it immediately follows that there exists the sign-Am^ degeneracy in measurement 
which combines any two of P, P T , and P s , and the explicit form of the degenerate solution 



can be obtained by the symmetry alone. Later in Sec. VII and Sec. VIII we will explicitly 
verify it by working out explicit solutions. 

A better understanding of the implication of the invariance may be achieved by drawing 
the bi-probability plot [21]. Here, we take the particular one in P — P T space [72] as shown in 
Fig. |2} The simultaneous invariance of P = P(^ M — > v^) and P T = P{y e — > v^) for any values 
of #13 and 5 means that one can find always two completely overlapping ellipses, one with 
positive and the other negative Am^. The 5 label, if placed onto the ellipses, are different 
between the two ellipses, and they are related by 5_ = tt — 5 + , where 5± denote the 5 for 
the respective hierarchies. More comments on the degeneracy in T-conjugate measurement 
will follow in Sec. IVII1 

Now, let us focus on the degeneracy with CP-conjugate measurement. Unfortunately, 
the similar simple-minded symmetry argument does not go through in the settings with CP 
and CPT conjugate measurement; The positive Am| x CP and CPT conjugate probabilities 



P CP and P CPT defined in Q and ([6j are not simultaneously invariant under (11). That 
is. P CP and P CPT are invariant under a transformation s — > a/ X + /X_s accompanied with 
the other transformations in (11). But, the transformation cannot be identical with (11) 
in matter. It is also true that the similar symmetry argument does not go through for the 
intrinsic and the #23 octant degeneracies. 
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FIG. 2: P-P T bi-probability plot. P - P T bi-probability plot by using the approximate formulas 
given in g and (||). sin 2 29 13 = 0.05 and NH (black solid), 0.061 and NH (red solid), 0.058 and 
IH (blue dotted), 0.069 and IH (blue dashed). 



However, we will show by the discussions throughout the following two sections and 



Sec. IX that there is a similar invariance of a pair of the oscillation probabilities, e.g., P and 



P CP in the CP-conjugate measurement, under the transformation 



si ->> s N = f N (si,5i,02™), 

5 1 5 N = r/ N (s!,5i,^ r 3 Ue ), 

Amjj — » +Am 31 , or - Am^, 

itrue , /itruc /ifalsc 



C -> or C- (12) 



In (12), si, <5i, and ^I;™ is the true parameters with subscript "1", the unique case with Ara- 
bic numerals. The degeneracy solutions are labeled by using Roman subscripts N= II— VIII. 2 
Alternative choices in the last two transformations determine whether the degeneracy is of 
the type involving the sign change of Am| 1; or the 823 octant flip. In this sense, the param- 
eter degeneracy is nothing but the statement of invariance of P and P CP under the discrete 



transformation (12). Furthermore, what is to be really stressed is that the mapping in (12) 



can be constructed by the basic three mappings. See Sec. [Vf 

III. INTRINSIC AND SIGN-Amjjj DEGENERACIES; CP-CONJUGATE MEA- 
SUREMENT 

A. Preliminary remarks 

Since this is the first section to actually solve the degeneracy problem to obtain the clone 
solutions we make some preliminary remarks. In the rest of this paper, we analyze the 



2 We note that each type of degeneracy is only two-fold merely because of the approximate form of the 
probabilities we use, or in other word, due to the smallness of #13. Existence of more solutions may be 
signaled e.g., by the sf 3 terms in the oscillation probabilities in a large-#i3 perturbation theory [73] . 
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structure of parameter degeneracy from various viewpoints, in particular, explicit analytic 
solutions, the symmetry aspect (as already mentioned), and use of the bi-probability plot to 
illuminate the respective characteristic features of the degeneracy. While the solutions were 
presented in a condensed way in [21] , we present a step-by-step derivation of the degeneracy 
solutions because it is reader friendly and makes the understanding of structure of degeneracy 
much easier. We emphasize that the analytic solution with the proper convention for its 
definition is the integral part of our discussion of solution network to be given in Sec. |Vj 
In most part of this paper, we confine ourselves into the oscillation channels between 
and v e and their antiparticles. In particular, we focus on v e (and z/ e ) appearance channel 
which will be available in conventional superbeam in this and the next sections. On the 
other hand, the T-conjugate channels, u e — > and v e — > z/^, would be provided by neutrino 
factory and the beta beam [7U [75]. Given understanding the degeneracy of the former 
channels, the corresponding informations on the latter may be obtained by regarding 5 by 



2n — 5. We include v T appearance channel in Sec. VIII 



In this paper we take the method for obtaining the degeneracy solutions for a given set 
of true parameters (s\, 5i). Equivalently, one can choose an alternative way of obtaining the 
degeneracy solutions as a function of "observable", e.g., P and P CP , as pursued in [29]. If 
one want to take this attitude, one can simply do it (at least numerically) by regarding that 
both the clone solutions and (P, P CP ) are parametrically represented by (si,5i). 

Notational comment: We denote the mass hierarchy of the true solution by the ± signs 
(+ for the normal and — for the inverted) to make the hierarchy choice always explicit. 
The relationship between the degeneracy solutions with input true mass hierarchies will be 
further discussed in Sec. IVl 

Now, let us start our discussion of parameter degeneracy by taking CP-conjugate mea- 
surement. The setting seems to be the most promising one experimentally in the near 
future. In this section we confine ourselves to the degeneracy solutions which have the same 
623 octant, though we treat generically the case of arbitrary values of 823 ■ The degeneracy 
solutions across octants of 6*23 will be discussed in the next section. 



B. The intrinsic degeneracy in CP-conjugate measurement 



With expression of the oscillation probabilities in ^ and (|4]), the intrinsic degeneracy 
solutions (si, 5i) (i=l, 2) in CP-conjugate measurement are defined with — > v e channel 
by 

P — Z = X±s\ + Y±si (cos S\ cos A 3 i =f sin 5\ sin A31) , 

P — Z = X±sl + Y±S2 (cos £2 cos A31 =1= sin 82 sin A31) , (13) 

and in CP-conjugate channel by 

P CP — Z = X T s± — Y T si (cos^i cos A 3 i ± sin^i sin A 3 i) , 
P CP — Z = X T s\ — Y T s 2 (cos 5 2 cos A 31 ± sin 5 2 sin A 31 ) . 



By subtracting two equations in (J 1 3|) and (|14|) respectively, we obtain 
Y ± 



Y+ 



s \ - s 2) + -TT- cos A 3 i(si cos^i - s 2 cos<5 2 ) =F ~rr~ srn A 3 i(si sin 8\ - S2sm8 2 ) 

x± x± 

Y Y 

— - cos A31 (si cos 5i — S2 cos #2) =F sin A31 (s% sin 5\ — s 2 sin #2) 

X T X T 



(si 



si) 



(14) 
0. 

0. (15) 
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From (15) we can obtain the expressions of cos <5 2 and sin <5 2 as 



where C^> is defined by 



2 

S 2 COS 82 = Si COS 5i ± 7 7) k (si — s?) , 

cosA 3 i(C(+)) 2 -(C(-)) 2 V 2 11 

2 C (+) 

s 2 sin 5 2 = si sin Si ± - — 7 ^ ^ (si — sf) , 

sinA3i(C(+)) 2 -(C(-)) 2 V 7 



c(±) = 2± ± 2L. 



(16) 



(17) 



- ^ S 2 — 1 gives a quartic equation for s 2 as (s 2 — s i)( s 2~ = 
••'•'" +: •• ' K " ••:♦•••■>+:-•'• "«» ; -— ^ to discussions of 



Inserting ( 16 ) into cos 2 <5 2 +sin 2 5 2 = 1 gives a quartic equation for s 2 as (S2 — s 
Of course, we obtain the trivial solution s 2 = si, the situation unique to d 
iegeneracy. The genuine intrinsic degeneracy solution is given by 



0. Of course, we obtain the trivial sc 
the intrinsic degeneracy. The genuine 

' 3 5i + Rsi sin 5i) (CH cos A 31 - RC {+) sin A 31 ) 



su = 



± 



II + R 2 ) 



(Si COS( 



+ 



1 



where R is defined by 3 



4(1 + R 2 ) 



(C*-> cos A31 - RC {+) sin A31) 1 



1/2 



(18) 



c(+) 

i? = ^py COt A31. 



Notice that the ± sign in (18) represents the mass hierarchy of the true solution. 
using (18) into (16) we obtain the solution of 8\\. From (16) one can obtain the followm; 
expressions: 



(19) 
By 



. 2Rsi I 

susm(5i + 5 u ) = T 

, - r> \ ( \ — R 2 \ 

sn cos (di + on) = 



'sin o"i + i? cos o\) (C H cos A 3 i - sin A 3 i) , 



Sl 



l + R\ 

=F on 7 D2 , (cos 81 - R sin 50 (C ( - } cos A 3 i - sin A 3 i) . (20) 

If we further expand the solution in (18) by the solar-atmospheric ratio ^"l 1 , we obtain 

sn = si + 2(1 ^ 2) (coso\ + sin 50 (C« cos A 31 - RC^ sin A 31 ) . (21) 

Similarly, the equation (20) also simplifies to cos (o\ + <5 2 ) = —(1 — i? 2 )/(l + R 2 ). These 
expressions reproduce the ones derived in [231 EH] • 



3 We note the relationship between the notations in this paper and in the previous papers 
here as MNP: C^lthis = 2CW| mnp, and i? I this = z 1 |mnp- 



denoted 
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C. The sign-Am 2 degeneracy in CP-conjugate measurement 



We turn to the flipped Am 2 -sign degeneracy in CP-conjugate measurement. The true 
input solution (si, Sx) and the opposite sign clone solution (s 3 , S 3 ) satisfy the following equa- 
tions. In the neutrino channel, 



P — Z = X±s\ + Y±si (cos 5i cos A 31 =f sin Sx sin A 31 ) 
P — Z = X T s\ + Y T ss (cos<5 3 cos A 3 i ± sin S3 sin A 3 i) 



(22) 



and in CP-conjugate channel 



P CP — Z = X T s 2 — F T si (cos^i cos A 3 i ± sin^i sin A 3 i) 
P CP — Z = X±s\ — Y±s 3 (cos<5 3 cos A 31 =1= sin<5 3 sin A 31 ) 



(23) 



By combining the first and the second equations in (22) and (23) we obtain 
X 4 



L ± „2 



X 3 



X 



s 3 + cosA 3 i 
=F sinA 31 



So — cosA 31 



— sicosdi - — s 3 cosd 3 



X 



Y- 



-Si sin 5\ + -r^Ss sin S : 



X: 



=p sin A 



31 



Yf c- Y± 

—s 1 cosb 1 - — s 3 cosd 3 

A± A ± 

—Si sin di + — s 3 sin dj 
A± A ± 



0. 



(24) 



Using defined in (17) we can simplify the equations. By subtracting and adding two 
equations in (24) we obtain 

T^ p - C (+) s 3 cos 63 cos A 3 i + C (_) s 3 sin 5 3 sin A 31 = 0, 



Tg£ -2s\ ± C ( ~ ) s 3 cos(5 3 cos A 3 i =f C (+) s 3 sin5 3 sinA 3 i = 0, 



(25) 



In (25), T±± and T-yf are defined as: 

Tgf = ±E (_) s? + D (+) sicos5icosA 3 i - D^Sisin^sin A 31 , 
T^f = E {+) s\ iD^sicosfJicosAai ^ £> (+) si sinc^ sin A 31 , 

where we have introduced the new notations and as 



2j(±) = 2± ± 

X_ X 4 



x_ x+ 



(26) 



(27) 



From ( 25 ) we obtain the expressions of cos 63 and sin S3 as 

1 



s 3 cos <5 3 
s 3 sin S 3 



1 



cosA 31 (C(+)) 2 - (<?(-))' 

1 1 
sinA 31 (C(+)) 2 - (C(-)) 2 
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[ C {+) t cp ±c H (T 2 c /-2s 2 )] 
[C^±CW (T 2 c /-2s 2 )] 



(28) 



We insert (28) into cos 2 £3 + sin 2 £3 = 1 we obtain the quartic equation for S3 as 

4 (1 + R 2 ) st - AU±sl + V± = 0, 

where 

21 



f/ ± = i cos 2 A 31 (C( 



± 



sin 



(29) 



V + 



cot 2 A 31 + 



c(+) 
c(-) 



(T : 



xT) 2 ± 



sm 



A 



31 



/ t cp t cp , / 1 



i? 2 ) (Tf/) 2 . (31) 



Notice that R is defined in (19). Equation (29) has the obvious solutions 



1 



2(l + i? 2 ) 



U±[±] 



Ul 



[l + R 2 )V ± 



(32) 



where [±]* denotes a temporary sign which is independent of the hierarchy sign. We discuss 
immediately below (Sec. HID) the way how to determine the sign convention. One can 
easily show by using Cauchy-Schwarz inequality that V± > 0, and hence the right- hand- side 



of (32) is non-negative. Therefore, there are four real solutions of S3 and the two positive 
ones are physical; The sign-Am^ degeneracy is two-fold. 

[I + R 2 )V± < defines the region in which there is 



The region specified by D 



sign 
± 



no sign-Am^ degeneracy solution. The region of no solution is displayed in Fig. |3jby taking 
the three typical setups for superbeam type experiments, SB1, MB1, and MB2 settings, 

4 In this figure the true mass hierarchy is taken to be the 



which will be defined in Sec. VIA 



inverted one. If we take the input normal hierarchy we must have the figure with 5 shifted 



by 7i, as one can confirm by looking at the white region in Figs. 10 and 12 in Sec. VI C This 



is expected by the general discussion to be given in Sec. |VB 

A notable feature is that the region of absent solution occupies mostly around 5 ~ tt/2 in 
the left two panels in Fig. [3] where L/E taken are at around the first oscillation maximum. 
The region is farthest to the "central region" populated by the both Am^-sign ellipses in 
the bi-probability space, and hence it is the region of lucky resolution of the sign-Am 2 ^ 
degeneracy |76j for the inverted (5 ~ 37r/2 for the normal) mass hierarchy. Generally 
speaking the no sign-degeneracy region grows for longer baseline, and the tendency continues 
to e.g., L = 4000 km and E = 20 GeV (Figs. PL8] and [19] in Sec. |VIE[). However, the feature 



changes for a region of the second oscillation maximum as seen in the third panel in Fig. |3j 
Because of the dynamic behavior of the bi-probability ellipses (see e.g., Fig. 2 in [53]) there 
are much better chance of having the sign-Am^ degeneracy. 



D. Problem of convention of labeling the degenerate solution 



We denote the two solutions in (32) as sni and srv- It is a highly nontrivial issue how to 



define these two solutions. In principle there are two ways: 



4 Apparently, a complete description of the no solution region of the sign-Am^ and the 823 octant degen- 
eracies seems to be lacking in the literature. 
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FIG. 3: Depicted as the shaded areas in the sin 2 2#i 3 — 5/tt space are the regions where no sign- 
Anijj degeneracy solution exists for the same three set of the baselines and the neutrino energies 
as in Figs. [T0| and |T2~1 The true mass hierarchy is taken to be the inverted one, and it may help to 
understand the relationship between the degeneracy solutions with different true mass hierarchies 
to be discussed in Sec. IV Bl 

Convention A: One can take the convention such that always Spy > Sm. That is, the plus 
and the minus signs in (32) correspond to Sjy and Sni, respectively. 

Convention B: One may choose the other convention such that the vacuum limit of the 
degenerate solutions can be taken smoothly. 

For reasons explained below we adopt the convention B. We note that D± sn denned as 
£,sign = u2_( 1 + R 2j y ± can be written as 



D s ± gn 



16 



cos 4 A 31 (C(-)) 4 {l-(^)' 



cot 2 A 31 (Tf/) 2 ± \ cot 2 A 31 C^C^T^ + \ cos 2 A 31 (^-)) 2 (l + R 2 ) T. 



CP 
2± 



.(33) 



In the vacuum oscillation limit, a ->• 0, CW, £> (+) , R, and all vanish, and #) = 2 
and DC") = CH = cU hold, where 



Ct) = lim 



2F+ 



2A 



21 



sin 2^i2 cot ^23- 



Then, D s l gn has a vacuum limit I)'"'' 



a-s>o X± sinA 3 i 

^oct-vac = Um ^ o ^sign _ ^ 

dT = C£> cos A 31 Qc?H cos a 31 ± Sl cos 5 



^ gn ) 2 where 



(34) 



(35) 



The smooth limit to the sign-Am^ degenerate solution in vacuum can be achieved by taking 
the sign convention 



6 III 



'IV 



1 



2(1 + ,R 2 ) 
1 

2(1 + R 2 ) 



D 



± 



(36) 
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In (36) we have taken the convention such that in a — > limit % and srv smoothly tend 



to s 



in 



and 



'IV > 



respectively, in vacuum defined in Sec. VII D Once the solutions of S3 



are specified with the well defined convention the solutions 6m and 5\y can be obtained by 



inserting s ra and s\y, respectively, into (28). 



We take the convention B because of number of desirable features. The matter pertur- 
bation theory [TTJ EE] can be formulated only with the convention because it requires the 
existence of smooth limit a — > in each solution. More importantly, the convention B makes 



the structure of the degenerate solutions transparent. That is, if we denote Sni in (36) and 



<5ni in (28) in an abstract fashion as 



sin 



then, one can show that 



srv 



mi 



/-CP sign 



6 



IV 



77± Ps1s >ia; 



CP sign/ r \ 



(37) 



(38) 



In this sense there is the one-to-one correspondence between the two intrinsic degeneracy 
solutions and the sign-Am^ degeneracy solutions, a charming property which one can enjoy 
only with the convention B. 5 The relations will be further extended into the other types pf 
degeneracies and completed in Sec. |V} 

However, there exists a somewhat disturbing feature of this convention; The solutions 
have discontinuity as a function of E, the neutrino energy. It by no means, however, that 
the solutions are unphysical. It merely implies that the two solutions interchange themselves 
at the discontinuous point. 



IV. 6 23 OCTANT DEGENERACY IN CP-CONJUGATE MEASUREMENT 

In this section we address the parameter degeneracy solutions, assuming that #23 7^ 7r/4, 
which have the different 823 octant from the true one. We denote the octant in which the 
quantity lives by the superscript "true" or "false" where the true 6*23 can be in either the 
first or the second octants. Our treatment of the octant degeneracy will be done under 
the approximation that the two solutions of 823 has the same value of sin 2823. In this 
approximation Y± is independent of the octant, but X± and Z in ([7]) have to have additional 
superscripts such as X± ne or X|* lse to indicate in which octant they live. Using the definition 
of X± and Z± given in (|7|) they are related with each other as 

X£ alse = cot 2 0* r 3 ue X£ ue , Z^ lse = tan 2 0* ue Z£ ue . (39) 



A. Intrinsic degeneracy in the false 623 octant 

We first discuss the case in which the intrinsic degeneracy solutions exist in a 823 octant 
different from the true one. The input solution (si, Si) and the different octant clone solution 



Notice that su > Si, or sn < Si, depending upon the region of experimental parameters. Therefore, once 



the correspondence relations (37) and (38) are established, we cannot take the convention A which implies 
that always Siv > Sm. If we try to solve the problem of obtaining the eightfold degeneracy solutions for 
a given set of (P,P CP ). it would be possible to take the convention A. 
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(s5, £5) satisfy the following equations 



P = X±s( + Y±Si (cos 61 cos A 3 i =f sin 61 sin A 3i ) + Z* 



£5 cos A31 =F sin 65 sin A 



31; 



7 false 



(40) 



in the neutrino channel, and 



P 
P 



;ytrue s 2 _ (cqs^ CO s A 3i ± sin ^ sin A 3 i) + Z t: 
J5s: feise s 2 _ ( cos ^ 5 cos a 31 ± sin <5 5 sin A 3 i) + Z f 



(41) 



in the CP-conjugate channel. The way we obtain S5 and S5 follows exactly the one for the 
sign-Am^ degeneracy solutions in Sec . Ill C Therefore, we can skip many equations and 
just say that (40) and (41) lead to 



s 5 cos 65 
S5 sin 5 5 



1 



cosA 31 (C<( +) )2 



sin A? 



31 ~ (Cf>) 



C^Tg P =F Ct ] (Tg P - 2sl) 



(42) 



In (|42j), we have used the new notations which essentially is a generalization of the previous 



one as well as new ones: 



G (±) 



false 
1 



J^false ' 
1 



ytrue 
F (±) = ±+ ± 



X 



false 



± 



false 



(2 



X false Jj£ false ' 
Hue ^false^ 



(43) 



and defined T£± and T£±, in parallel with and in (26), as 



nCP 

■3± 



±F { - ) s\ ± G H + Ci +) S1 cos 5i cos A31 - C { f ] Sl sin ^ sin A 31 



(-), 



T£ p = F (+) s 2 + ± Ci -) si cos5i cos A 31 =f c! 2 +} s 1 sin Si sinA 31 



(+), 



Inserting (42) into cos 2 5<s + sin 2 5 5 = 1 we obtain the quartic equation for s 5 as 

4(1 + R 2 2 ) sj - m± 8 \ + J± = 



(44) 



(45) 



which is actually a quadratic equation of s 2 3 because of the quadratic dependence on s 13 
of S5COS55 and S5sin5 5 . Thus, there exist only two physical (i.e., positive) solutions, which 
implies that the octant degeneracy is two-fold. In (45) H and / are defined as 

2 



H± = -cos 2 A 31 (C 2 



("h2 



1 



a 



(+) 



a 



(-) 



sin 2 A 31 [c^ 



T 3± + (1 + Rl) Tfl, (46) 



cot 2 A 



31 



C. 



(+)' 
P) 



( T 3±) 



sm 



2 

^A 



31 



C, 



(+)' 
H 



^Tif + (1 + i? 2 ) (Tf/) 



! .(47) 
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where R2 = -f^ycot A31. Then, the octant degeneracy solution s 5 is given by 



(+) 



2(1 + RV, 



H±[±] 



oct-intr 



(48) 



where [±]* is the temporary sign to be specified below, and £)° ct - mtr is defined as 



^oct-intr 



a 



(+)' 



16 



cos 4 A 31 (C 



HN4 




cot 2 A 31 (T, 



cot 2 A 3 id +) C. ( - ) 



T, 



CP 
3± 



cos 2 A 31 (cH) 2 (i + it:i)rf/ 



(49) 

The region defined by £)° ct - mtr < defines the region in which there is no intrinsic degeneracy 
solution in an octant different from the true # 23 . Once there is a solution, it must be obvious 



that s 2 is positive definite following the similar argument as in Sec. Ill C 



The region of no #23 octant degeneracy solution is displayed as the white regions in the top 



two panels in Figs. 14 and 16 for superbeam type settings, SB1, MB1, and MB2 (to be defined 
in Sec. VIA). The corresponding informations for neutrino factory setting NF are given in 



Figs. 18 and 19 In fact, one can observe that the feature of no octant degeneracy region is 
insensitive to the baseline and energies, but depend on in which 6* 23 octant the true solution 
exists. If #2 3 ue is i n the first octant, no-degeneracy regions are around 5 ~ tx/2 and 37r/2, 
whereas if it is in the second octant it is confined into the small #13 region, sin 2 2# 13 < 10~ 3 . 
Both of the features can be easily understood by drawing the bi-probability plot; The region 
of no degeneracy solution is the one spanned only by the ellipses with a single octant #23 

us. 

Now, we have to revisit the issue of convention to define unambiguously the octant 
degeneracy solutions. We take the following new convention: 



Convention C: We define (sy,5y) and (syi,5vi) such that they have a smooth limit to 
the intrinsic degeneracy solutions (si,5i) and (sn,5n), respectively, when the maximum #23 
limit #23 — > it /A is taken. It can be understood as a consistency condition. 

The convention allows us to formulate perturbative framework with use of the small expan- 
sion parameter #23 — 7r /4, as we will do in Appendix [Bj Notice that the degenerate solutions 
(siij^ii); (sin, 5m), and (spy, <W) themselves obtained in Sec. Ill are valid independent of 
the value of #23. 

In the maximum 6 23 limit, FH and all vanish, and F« = 2, C { 2 ±] = and 

R2 



R hold. Then, £)° ct - mtr has the maximum 623 limit 



£)Oct-intr-max __ Itjq 

6»23^-7r/4 



£)Oct-intr 



= (d° Ct " intr ) 2 



(50) 



where 



^oct-intr 



cosA 3 i < 1 



\c(-)J 



x 



cos A 31 (C ( 



Ci 1 



± C*- ' (si cos^i + Rsis'm.81) 



(51) 
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The smooth limit to the same-octant intrinsic degeneracy solution can be achieved by taking 
the sign convention 



'v 



6 VI 



2(1 + R 2 ) 



2(1 + Rl) 



H. 



± 



^oct-intr 



£)oct-intr 
^oct-intr)2 



H± + d% 



'oct-intr 



^oct-intr 
^oct-intr ^ 2 



(52) 



One can easily verify that in the maximum 823 limit sy and Syi tend to s\ and su, respectively 



B. Sign-Am 31 degeneracy in the false #23 octant 



We next discuss the sign-Am^ degeneracy in a 823 octant different from the true one. 
The input solution and the false octant clone solution (s 7 , 5 7 ) satisfy the following 

equations 



P = X±sf + Y±si (cos 8x cos A 3 i =F sin 5 X sin A 31 ) + Z 
P = X^ lsc s 2 7 + Y T s 7 (cos 5 7 cos A 3 i ± sin 8 7 sin A 31 ) + Z 

in the neutrino channel, and 



true 



false 



(53) 



p cp = x^sj - Y T si (cos5i cos A 31 ± sin^ sin A 31 ) + Z t 
P CP = X^ lse s 2 - Y±s 7 (cos 5 7 cos A 31 t sin 5 7 sin A 31 ) + Z 



false 



(54) 



in the CP-conjugate channel. Proceeding along the same way as in Sees. Ill C and IV A| we 
obtain 



s 7 cos S 7 
s 7 sin S 7 
where we have defined 



1 



1 



cosA 31 (c<+))2 _ (C { f y ) 2 



sin A? 



31 (& 2 +) y - (cf) 



Cfhtl p ± Ct ] (iff - 2s 2 ) 

df>Tg±dp (iff -24 



(55) 



T°£ = ±E i f ) si t + D ( 2 +) Sl cos 6^0$ A 31 - D { 2 ^ sint^ sin A 



31) 



r CP 
■6± 



E { 2 +) s\ + ± D { ~ ] Sl cos^i cos A 3 i =f D { 2 +) s x sinc^i sin A 31 . 



i(-) 



(+), 



(56) 



We have used the notations 



D 



(±) 



G <±) 



J£ false 
1 



J£ false ' 
1 



E. 



(±) 



X 



true 



± 



X 



true 



± 



J£ false J£ false ' 
1 rue ^false^ 



J£ false J£ false 

Inserting (55) into cos 2 5 7 + sin 2 5 7 = 1 we obtain the quartic equation of s 

4 (1 + R 2 2 ) s 7 - U±s 2 + K± = 0, 



13, 



(57) 



(58) 
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which is actually a quadratic equation of s 2 3 because of the quadratic dependence on s i3 
of S7COS57 and Sysin^. Thus, there exist only two physical (i.e., positive) solutions, which 
implies that the octant degeneracy is two-fold. In (58) J and K are defined as 



cos 2 A 31 (Cf >) 2 





sin 2 A 



T 5± + (1 + Rl) T 6 C ± P , (59) 



cot 2 A31 + 



a 



a 



(-) 



a 



(+)' 



sm 



V31 



Then, the sign-Am 31 degeneracy solution across 823 octant is given by 

1 



■4[±] 



^oct-sign 



7 2(l + i? 2 ) 

where [±]* is the temporary sign to be specified below and £)° ct ~ sl § n j s defined as 



(61) 



^oct-sign 




1 

16 



cos 4 A 31 (C, 



(-h4 




cot 2 A 31 (Tf/) 2 ± \ cot 2 A 31 C«C«Tf/ + I cos 2 A 31 (C^) 2 (l + i? 2 ) T ( 



CP 
6± 



(62) 

The region defined by £)° ct_sl s n < q defines the region of no sign-Am 31 degeneracy solutions 
in the false #23 octant. The region of no sign-octant degeneracy solution is displayed as the 



white regions in the bottom two panels in Figs. [14] and [16] for superbeam type settings, SB1, 
MB1, and MB2. The corresponding informations for neutrino factory setting NF are given 



in Figs. 18 and 19 Once there is a solution, the same argument as before assures that the 



solutions for s 2 in (|61|) are positive definite. 



To define unambiguously the sign-Am^ degeneracy solutions in the false #23 octant we 
need the following new convention. That is, we need to take both the maximum #23 and the 
vacuum limits. 



Convention D: We take the convention such that (svn, #vn) and (syin, 5viii) tend to 
(siff, Sjjf) an d (sjv°, <5Jv c )' respectively, in the simultaneous maximum-^ and the vacuum 
limit. 



One can easily show in the combined limit that limg 23 . 



^oct-sign 



^(023 



7r/4, vacuum 



^oct-sign = ( rf oct-sigu )2) where 

± 1^23 — tt/4) with d± gn being defined in (35). The smooth limit to the maximum 
# 23 of the sign-Am 31 degeneracy solutions can be achieved by taking the sign convention D. 
The solutions read 



'VII 



,2 

6 VIII 



2 (1 + Rl) 



2(1 + m) 



J4 



^oct-sign 



^oct-sign 
(d° Ct - Sign ) 2 



t 1 joct-sign 



^oct-sign 
(4 Ct " Sign ) 2 



(63) 
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V. STRUCTURE OF PARAMETER DEGENERACY 



We now make the structure of parameter degeneracy transparent based on knowledges 
obtained in the previous two sections by using CP-conjugate measurement. First, we sum- 
marize the relationships between the true and the degeneracy solutions. 



A. Mappings between the true and the degeneracy solutions 



Let us start by putting the relationship between each intrinsic degeneracy pair of the 
solutions in order. If we denote the relationship between the intrinsic degeneracy solution 



derived in Sec. Ill B as sn = C±( s i> #23"°) an d = ^±( s i, ^23 e ) then the other intrinsic 
degeneracy pairs of the solutions satisfy 



S PV — C=f( s IH' ^IIIj #23 Ue ); 

svm = Ct( s vii, "vii, 7r/2 - 6>23 UC ) 



"TV = <-=p(SlH, Sm, ^23°) 

Syi = t±(s V , Sy, 7r/2 - 6» 23 
$vm = L T (syii,5yii,7r/2 



truc\ 



/itrueN 
^23 J> 



(64) 



completing the one-to-one correspondence between them. We then summarize the one-to- 
one correspondence relations between the solutions which involve the Am^-sign and/or the 
#23 octant flips: 



sin = 

sy 



_ c CP oct 



si,8x,8 23 j 



syii 



,-CP oct-sign/ r n 
4± {Si,di,V 23 



SVII 



/-CP sign 



true 
23 
true 



■SVII — s T 



Sy,5y,Tl/2 - 

tCPoct/ e c 
C,m [Sill, 



/itrue ' 
,0111, &23 . 



Sy = i 



CP oct-sign/ r ntrue- 



SlV = 
SVT = 
SVIII 
SVIII 
SVIII 

svi = 



/-CP sign/ c ntrue\ 

?± lsii,0ii,y 23 ) 



^CP oct 1 



SlI, $11, #23 



5^23 11 
true\ 



CP oct-sign/ r ntrnc\ 

{sn, On, V 23 ) 



/■CP sign/ r 
f± {Syi,8 

^CP oct 1 



VI 



tt/2 



6T) 



s Wi 8iy, 9. 



/■CP oct-sign/ r 

i T [sw,o 



23 
IV' # 



true\ 
> "23 ) 



(65) 



The functional form of £± P Slgn , £± p oct , 



and £ ± P oct_slgn are defined in (|37|) in Sec. 



first line in (52) in Sec. IV A, and the first line in (63), respectively. There exist t 



HID 



the 
re similar 

relationships between <5's through the function rf^ Slgn , but we do not display them explicitly 



here. The correspondence relations (65) between the true and the degeneracy solutions 



II— VIII are pictorially represented in Fig. HI The relationship between the degeneracy 



solutions in (64) and (65) as well as in Fig. \4\ is the precise meaning of the statement that 



nature of the parameter degeneracy is the intrinsic degeneracy duplicated respectively by 
the Amnj-sign and the #23 octant flips. 



It is easy to prove (|65|) by considering the original defining equations for the degeneracy. 

s x and tj2 P Slgn (s v , Sy, tt/2- 



For example, it is easy to show that Slgn (sv, Sy, tt/2 — 6. 
5x are solutions of the equation 



/itruc 
U 23 



true "| 
23 ) 



Xf 1Se) 4 + Y ± Sy C0S(CV ± A31) + ^ (falSC) 
X (false) s 2 _ y ^ Sy cos(Jv t Asi) + ^(false) 



X^ lse) si + Y T s x cos(5 x TA 31 ) 



^ (false) 



X^s 2 x - Y ± s x cos(5 x ± A 31 ) + 



(66) 



Note that left hand side of Eq.(66) can be expressed by s\, Si, and #23™ from Eq.(|40j) and 
(41), it is nothing but the equations of sign-octant degeneracy. 
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FIG. 4: The correspondence relations (65) between the solutions, the true one and the clones 
II— VIII, are pictorially represented. 



Though we do not present explicit formulas the similar structure exists in all the degen- 
eracy solutions in other settings, T-conjugate, Golden-Silver, and CPT-conjugate measure- 
ment to be discussed in the following sections. 

B. Relation between cases of true normal vs. true inverted mass hierarchies 

Here, we note an important property of the degeneracy solutions. Namely, if we know 
the degeneracy solutions for the input true normal mass hierarchy, then the solutions for 
the true inverted mass hierarchy can be immediately obtained from the former ones. 6 Since 
the statement can be confusing to the readers, we give an explicit proof as below. 

Let us take the sign-Am^ degeneracy (without 9 23 octant flip) for definiteness. For 
clarity, we denote the degeneracy solution for the case of true normal and true inverted 
mass hierarchies as (s 3 n,5 3 n) an d (s 3 i,5 3 i), respectively, where the subscript "3" implies 
either III or IV. For a given set of the probabilities P and P CP , assuming that the true mass 
hierarchy is normal, the true solution (si,5i) and the fake one (s 3 n,#3n) satisfy 

(P =) X + s\ + Y + s x cos(8 x + A 3 i) + Z = X_sl N + F_s 3N cos(5 3N - A 31 ) + Z, 
(P CP =) X_s\ - Y_ Sl cos^ - A 31 ) + Z = X + s 2 3N - Y + s 3N cos(5 3N + A 31 ) + Z. (67) 

If the true mass hierarchy is inverted, then the degeneracy solution satisfies a different set 
of equations as 

(P =) X^s\ + cos(5i - A 31 ) + Z = X+4 + Y + s 31 cos(5 3 i + A 31 ) + Z, 
(P CP =) X + s\ - Y +Sl cos^ + A 31 ) + Z = X_4 - F_s 3I cos(5 3I - A 31 ) + Z. (68) 

We define <5i ne w and <5 3 i new as 5\ = <5i new + vr and <5 3 i = S 3 i new + ir, respectively, and rewrite 



(68) by using them. It reads 

Xs\ - cos(5i ncw - A 3 i) + Z = X+s 3I - Y + s 3 i cos(<5 3Tncw + A 3i ) + Z, 

X + s\ + Y + s 1 cos(5i new + A 3 i) + Z = X_s 3I + F_s 3I cos(5 3 i ncw - A 31 ) + Z. (69) 

6 In spite of this charming and useful property we have been denoting explicitly the input mass hierarchies 
as the ± signs in the subscript for clarity and simplicity of keep tracking of the true mass hierarchy we 
are working. 
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Comparison between (67) and (69) tells us that if the set (s3n,<53n) is the solution to the 
sign-Am^ degeneracy equation for the true normal hierarchy, then the set (s 3 i, ^inew) is the 
solution for the true inverted hierarchy, provided that the true value of 5(— 5i) is replaced 
by 5inew Stated more explicitly, if the first two equations in (12) define the sign-Am^ 
degeneracy solution for the true normal hierarchy, 



(70) 



then, the degeneracy solution for the true inverted hierarchy is given, if expressed in terms 
of the true input parameters, as 



«3i = fafai, Si-n), 6 3 j = r] 3 (si, 6t - 7r) + n. 



(71) 



In other word, the mapping functions £± and r)± are related with each other as 



7T 



7T. 



(72) 



Thus, the sign-Am^ degeneracy solutions for the true inverted hierarchy are essentially 
determined by the solutions for the true normal hierarchy. The relationship between the 
degeneracy solutions with the normal and the inverted hierarchies as the true solution may 
be understood better by comparing the no-slotuion regions given in Fig. [3] which are drawn 



with the inverted hierarchy as the input true solution to those of Figs. [UN and 12 



It can be easily seen that the same treatment goes through for other types of degenera- 
cies, the intrinsic and the sign-Am^ degeneracies in the same or the different 623 octants. 



Therefore, the relationship (72) holds for the whole eightfold degeneracy. 



C. Asymptotic expansion 

At the end of this section, which is devoted to illuminate general properties of the de- 
generacy solutions, we make comments on high-energy behavior of the solutions. It would 
help clarifying some features of the energy dependence of the degeneracy solutions which 
are discussed in the next section. 

First of all, it can be easily verified that all the degeneracy solutions sn and <5n 



(N=II— VIII) derived in Sees. Ill and IV have finite asymptotic limit as E — > 00. Then, we 
note an interesting property that they are invariant under the transformation E — > —E, or 
A ji —> — Aji (j = 2,3). The invariance, of course, stems from the one of the oscillation 
probabilities. Though the transformation is unphysical in nature, the transformation prop- 
erty is useful. It means that when we do asymptotic expansion of the degeneracy solutions 
as sn = S n= o a n (^3i) n the odd terms are absent: 

SN = + a^Aj, + O (A*J , 

6 N = b^ + b^Al 1 + (Aij. (73) 

Absence of the first-order term in 1/E implies that onset to the high-energy asymptotic 
behavior of the degeneracy solutions is relatively fast, as we will confirm in the next section. 
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VI. OVERVIEW OF THE EIGHTFOLD PARAMETER DEGENERACY IN CP- 
CONJUGATE MEASUREMENT 



In this section, we try to give an overview of the intrinsic, the sign-Am^, and the #23 
octant degeneracies. In fact, the features of the degeneracy solutions are quite different 
for differing baselines and neutrino energies, which makes the overview in a genuine sense 
extremely difficult. Therefore, we restrict ourselves in this paper into a few typical settings 
which may be relevant for the settings of future neutrino experiments discussed in the 
literatures. If the readers want to examine features of the degeneracy with some alternative 
experimental parameters, they can do it quite easily by using the analytic solutions presented 
in this paper. 

We want to warn the readers that all of our comments to be made in this section are 
qualitative in nature. Therefore, when we say, "spectrum analysis would resolve the degen- 
eracy" it actually means that it may be possible to resolve it if appropriate experimental 
settings are provided. 7 Yet, we try to be based the experiences gained in some previous 
analyses. The readers may still wonder whether the discussion of degeneracy based on the 
probability makes sense because the observable in the experiments must be obtained after 
convolution with neutrino fluxes and cross sections. However, this is not the only possible 
attitude to take. One can, in principe, obtain the "experimental data of probability" by 
de-convoluting the fluxes and cross sections as shown in [SUj . 

We also remind the readers of the fact that limited statistics in measurement has a non- 
trivial impact on the features of the degeneracy to be observed in the actual experiments. If 
two (or more) degeneracy solutions exist in nearby locations in the mixing parameter space 
they can merge together, producing an apparent single solution, which could be misinter- 



preted as no degeneracy. We have discussed in Sees. [ITT] and [TV] the regions of no degenerate 
solutions. However, there might be cases that the degeneracy solutions do exist in these 
regions because of lack of statistics and/or shift of the allowed regions due to systematic 
uncertainties. Discussion of these features is outside the scope of this paper. 



A. Variables used for display and baselines and neutrino energies adopted 

We try to illuminate some characteristic features of the degeneracy by presenting the 
differences between the true solution and the fake ones. To display the difference between 
the solutions we define the ratio as 



sin 2 26^ - sin 2 26\T 

itrii' 
; 13 



sin 2 29 mc 



where N = II- VIII denote the degeneracy solution labels. We use the variable sin 2 29^ 
because it is closer to the experimentally measured quantity. Similarly, we define the quantity 
D N to represent the differences between the true and the clone solutions. For this purpose, 

7 Similarly, when we say "the degeneracy A is easier to lift than the degeneracy B" it actually means so 
provided that an appropriate experimental condition is prepared such that the similar sensitivities would 
be expected for both the solutions A and B. 
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there are two appropriate ways to define it, the types (1) and (2); 



4 X) - { — \ 4 2) - - ~ { 6 ] - (75) 

7T 7T 

We use either one of or whichever appropriate depending upon the degeneracy 
types. 

As typical experimental settings, we use the following four cases of baselines and neutrino 
energies: 

• SB1: L = 300 km, E = 700 MeV; A short baseline low energy (and v^) superbeam 
near the first oscillation maximum 

• MB1: L = 1000 km, E = 2.33 GeV; A medium baseline superbeam near the first 
oscillation maximum 

• MB2: L = 1000 km, E = 780 MeV; A medium baseline superbeam near the second 
oscillation maximum 

• NF: L = 4000 km, E = 20 GeV; A typical setting for neutrino factory 

The first three settings, SB1, MB2, and MB1, are examined in the following three subsections 



(Sees. VI B VIC, and VI D), while the last one, NF, in Sec. VIE We should remark here 
that the energies corresponding to each baseline are chosen rather arbitrarily just for display, 
while their order of magnitudes are dictated by the baseline distances. 8 In some limited cases, 
the features of degeneracy solutions of the first two cases, SB1 and MB1, are so similar that 
we omit MB1 plots. Comparison between MB1 and MB2 settings, former (latter) being 
around the first (second) oscillation maximum, would be interesting to know physics behind 
the potential of the two-detector setting [81] and/or the BNL-type wide band beam approach 
[52]. 

The key to resolve the degeneracy is to utilize spectrum informations. Therefore, we 
also present the energy dependence of the difference between the true and the degeneracy 
solutions. A point of interest is how the energy dependence differs among the three different 
types of the degeneracies. 

We note that in all the figures presented in this section we take the normal mass hierarchy 
as the input true solution. If one wants to have the corresponding informations for the 
inverted mass hierarchy, one can do it just by changing the ordinate label of the figures as 
5 true — 7r not <5 tme . This is discussed in detail in Sec. VB 



We use the following values for the mixing parameters as summarized below: Am^ = 
7.9 x 10~ 5 eV 2 , sin 2 6» 12 = 0.31, and Ami = 2.5 x 10~ 3 eV 2 . The matter density is taken as 



p = 2.8g/cm 3 for SB1, MB1, and MB2 settings, and p = S.Gg/cni 3 for NF setting. 



To avoid the energy of perfect oscillation maximum of about 600 MeV at 300 km, where the features 
of the degeneracy are special, we tentatively added 100 MeV. Hence, e, the deviation of A31 from ir/2 
is given by A 31 (L = 300 km , E = 700 MeV) = tt/2 - e from which e is determined as e = \\ = 0.224. 
For the second oscillation maximum, we have arbitrarily chosen that A3i(L = 1000 km , E) = 3(ir/2 — e), 
that is E = 780 MeV. 
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B. Intrinsic degeneracy in the true #23 octant 



In Figs. 5 and 6 we present Ru, the normalized difference of sin 2 2# 13 , and D 



(2) 
II > 



difference of S/n defined in (75), respectively, between the true and the intrinsic degeneracy 
solutions for the two typical cases of energies and baselines, SB1 and MB2. We do not 
present the same plots for MB1 setting because they are so similar to those of SB1. A color 

(2) 

variation is used to clearly represent the ratio Ru and in a visual way, which will be 



used later also to all of i?N and D$. From blue to red Ru and D\^' vary from —1 to +1. 
The only exception to this rule is R^ at color graduation of the deepest red; It contains the 
region with R^ greater than 1. Notice that there is no region of i?N < — 1 by definition in 
(74). 9 We should note that in the case of D$, unlike R^, the deep blue region smoothly 



)(2) 



continues to the deep red because of the periodicity in S. 
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FIG. 5: The ratio R u = [sin 2 26^-sin 2 29\f c ] / sin 2 29\f e defined in (74} is presented in sin 2 29^ c - 
5 tvue /ir space for the two typical cases of energies and baselines, SB1 and MB2, defined in Sec. VIA 



We do not present the same plots for the setting MB1 because they are so similar to those of SB1. 
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FIG. 6: The normalized difference = [5 U — (tt — <5 true )]/7r, defined in (75) is presented in 



sin 2 20f 3 ue 

not present the same plots for the setting MB1 because they are so similar to those of SB1. 



5 ue /7T space for the two typical settings SB1 and MB2 defined in Sec. VIA We do 



9 Here we need to mention about some details of color variation. 20 color graduation are used to draw i?N 
and D$ so that a single color graduation spans 5% of the entire region. 
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FIG. 7: Bi-probability plot for settings of SB1 (left panel) and MB2 (right panel). In both plots 
the true value of #13 is taken to be sin 2 2^13 = 10 -3 . The ellipses of the intrinsic degeneracy 
solutions are depicted by the red solid and red dashed lines. 



One of the most notable features in Fig. [5] is a clear difference between SB1 (left panel) 
and MB2 (right panel) settings. In large 6*13 region in SB1 setting, sin 2 26*i 3 > 1CT 2 , Ru is 
small. At small #13 in SB1 setting and at large #13 in MB2 setting, Ru is large and positive 
(negative) at 5 ~ (7r). In small 13 region in MB2 setting, sin 2 26*i 3 < 10~ 2 , Ru is large 
and positive independent of 5. 10 Let us understand these features. 



We start from the above first feature. It can be understood by the analytic solution (18) 
If #13 is relatively large, S13 ^> Am^/Am 2 ^, s\ is the dominant term in su in (21), and 



Ru ~ O (sT^Am^/Am 2 ^) . Therefore, Ru is small at large 6*13 in SB1 setting. Now, the 
behavior of Ru at small #13 can be easily understood by looking into the bi-probability plot, 
the left panel in Fig. [7j The degeneracy ellipse which shares the point around 5 = (5 = ix) 
of the true ellipse is the dashed (solid) one with considerably larger (smaller) 6* 13 . The similar 
consideration explains the feature of Ru at large #13 in MB2 setting. The remaining feature 
that needs explanation is the large positive Ru at small # 13 in MB2 setting. At such small #13 
as sin 2 2#i3 ~ 1CT 3 and the baseline L = 1000 km, the oscillation probability is dominated 
by the solar term Z. Since it is independent of S the probability ellipse shrinks to a small 
"circle" , as can be seen in the right panel in Fig. [7} Then, the degeneracy solution ellipses 
are inevitably large as indicated by the red solid and dashed lines, resulting degenerate 
solutions much larger than the true S13. 

We observe for D$ plotted in Fig. [fi] that in SB1 (and MB1) setting Djp is small in 
the entire region of sin 2 2# 13 — 5/n space covered. 11 It means that the approximate formula 
8 ~ 7T — <5 true works well |23j . It is a nice feature of measurement by SB1 setting because 
CP violation is unlikely to be confused with CP conservation. Whereas for MB2 setting 
Djj is small only in a limited region —0.2 < 5/n ^ 0.2, and in a small strip around 5 ~ 
at large #13. The deviation from the approximation 8 11 ~ 71 — <5 tme is significant in the 



10 



At extremely small #13, the left plot in Fig. [5] for SB1 almost looks like the right plot for MB2, but with 
scale of sin 2 2£>i3 two orders of magnitude smaller than that of the right panel in Fig. [HJ 
11 Hereafter, when we talk about sin 2 2^13 and S such as sin 2 2#i3 — 8 /it space, it actually means the 
sin 2 26**3 Ue and S true , respectively. We use the simpler notation to avoid cumbersome superscript "true" as 
much as possible. 
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second and the third quadrants of S, in particular in region sin 2 2#i3 < a few x 10 2 in MB2 

setting. It is possible to understand this behavior of D^' qualitatively at very small #13, 
sin 2 26 l 13 ~ 1CT 3 by using the bi-probability plot for MB2 setting (right panel in Fig. For 
the true value of 5 truc ~ the degeneracy ellipse is depicted by the solid line and 5n o± n, 

( 2") 

which implies <C 1. On the other hand, for 5 tmc ~ n the degeneracy ellipse depicted 

(2) 

by the dashed line touches to the true ellipse also at around 5 ~ tt, hence D\ Y is of order 
unity. Notice again that the deep blue region smoothly continues to the deep red because 
of the periodicity in S. 
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FIG. 8: The energy dependence of sin 2 2#}g is plotted for the two typical settings SB1 (left panel) 

The true value of #13 is taken as sin 2 2#i3 = 0.05, 
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and MB2 (right panel) defined in Sec. 
which is indicated by the horizontal solid line in the figure. 
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is plotted for the two typical 



(2) 

FIG. 9: The energy dependence of the ratio defined in 
settings SB1 (left panel) and MB2 (right panel) defined in Sec. VIA The true value of #13 is taken 
as sin 2 26 13 = 0.05. 



One of the most important issues to find possible ways to resolve the degeneracy is to 
know the energy dependence of the difference between the true and the clone solutions. 
Therefore, we present in Fig. [8] the energy dependence of sin 2 29^ for SB1 (left panel) and 
MB2 (right panel) settings, assuming the true value of sin 2 26*13 = 0.05. It may be regarded 
as a typical value for relatively large #13 to which we will have an access by the ongoing 
experiments. Hereafter, whenever we present the similar figures of the energy dependence 
of the degeneracy solutions, we use only the three values of 5, 5 = 0,n/3, and 27r/3, for 
better visibility. (We have tried S = 0,27r/3, and 47r/3, but the latter two curves tend to 
overlap.) Similarly, in Fig. M the energy dependence of is plotted for the same settings, 
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SB1 and MB2, with the same true value of #13. Again, we do not present the case of MB1 
setting because the plots are very similar to those of SB1 apart from minor differences at 
low energies, E < 1 GeV. 12 

Clearly, there exists a significant energy dependence of sin 2 26^ even for SB1 setting. 
One can see that they vary by a factor of 2—4 (30%— 40%), or more at low (high) energies 
depending upon 5 in region of E = 0.4 — 1.2 GeV for L = 300 km. It must be contrasted 



to almost flat curves of energy dependence given in the following subsections, Figs. [TT 
(Sec. VIC) and 15 (Sec. VI D) for the sign-Am^ and the octant degeneracies, respectively. 
Then, the spectrum analysis must be powerful in resolving the intrinsic degeneracy. It has 
been seen to be the case in the analysis of T2K II experiment [13] done in [53]. For D^' the 
energy dependence is significant only at low energies, below the first oscillation maximum, 
where usually the signal-to-background ratio is not helpful. Therefore, potential power for 
the spectrum analysis relies more on the energy dependence of sin 2 20*3 not on 5ns. 

The energy dependences of sin 2 26^ 3 and plotted for MB2 setting magnify the low 
energy part of MB1 setting (not shown) at the same baseline of L = 1000 km. The energy 
dependences are far more pronounced and depend sensitivity on 5, and at some particular 
energies sin 26f 3 and are pinned to sin 2 26\f e and zero, respectively. To understand 
better these features, we first note that the oscillation probabilities show violent energy de- 
pendences below the first oscillation maximum. The (^-dependent strong energy dependence 
of the degeneracy solutions, however, is under restriction by pinning to the true value or 
zero at the energies corresponding to the n-th oscillation maxima (n = 1,2,3), as seen in 
the left and the right panels of Fig. [8j In SB1 setting (300 km) the first oscillation max- 
imum is reached at E ~ 600 MeV, while for L = 1000 km the first, second and the third 
minima are approximately at E ~ 2 GeV, 700 MeV, and 400 MeV, respectively. The zeros 
of sin 26f 3 - sin 2 26\^ e and arise due to the special feature of the parameter degener- 
acy at the n-th oscillation maxima [271 EE] , which can be seen explicitly from the formulas 
in Sec. HI B It can also be intuitively understood by having a flattened ellipse in the bi- 
probaility diagram at the oscillation maxima [24j. 13 The regularity of alternating "pinning 
to zero" and violent energy dependence gives us at least clear picture of what is seen in 
energy dependence of the degeneracy solutions in MB2 setting. The similar features will 
be seen in many figures of energy dependence in MB2 setting presented in the rest of this 
paper, which will allow analogous explanations. 

Despite a possibility of confusion due to too complicated dependence on energy and 5 
at MB2 only setting, combination of MB2 at somewhat off the oscillation maxima with 
more quiet SB1 settings would be an ideal machinery for resolving the degeneracy. This 
was observed to occur in [531 El] which utilizes the informations at the second oscillation 
maximum by a Korean detector, and most probably gives an explanation for high sensitivity 
achievable in the BNL-type wide band beam strategy [521 157] . 



12 For more extensive presentation of these plots including those of MB1 setting and energy dependence 
plots with other values of S, see [75], and partly |82j . 

13 The similar zero at E = 1 GeV can be understood as shrinking the ellipse into a small size at the 
oscillation minimum and is not interesting to us. 
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C. Sign-Am^ degeneracy in the true #23 octant 



We now turn to the sign-Am^ degeneracy which exists in the same #23 octant as the true 
one. Since there are two solutions, (%, 5m) and (srv, 5iy), we present them in the same 
figures. In Fig. 10, Rm and Riy defined in (74) are plotted in sin 2 26> 13 — 5/ir space. White 



region is the region of no degenerate solution as discussed in Sec. |III C[ We now show also 
MB1 case because difference from SB1 setting becomes non-negligible, in particular in small 
9 13 region. Nevertheless, the difference is not so significant at large #13 sin 2 2#i 3 > 1CT 2 
apart from the change in the no-solution region. 
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FIG. 10: The ratios R m 



[sin 2 26™ 

three panels) defined in (74) in sin 2 26 



sin 



13 



]/ sin 2 

5/tt space is presented for the three typical settings SB1 



26» t 1 r 3 ucl 



28\™ e (upper three panels) and Riy (lower 



(left panel), MB1 (middle panel), and MB2 (right panel) defined in Sec. VIA The regions of white 
color denote the regions of no sign-degeneracy solution. 



By comparing Fig. 10 with Fig. [5J it is evident that the difference of sin 2 2#i3 between 
the true solution s\ and % is much smaller than the case of intrinsic degeneracy solutions 
for SB1 and MB1 settings. Furthermore, the energy dependence of sin 2 2#} I 3 I is much milder 
than the case of solution (s n , 5u) of the intrinsic degeneracy, as one can clearly see by 
comparing Fig. 1 1 with Fig. [81 These features make resolution of the sign- Am^ degeneracy 



much more difficult compared to the intrinsic degeneracy in these settings. 

The difference between energy dependences of the intrinsic and the sign-Am^ degen- 
eracies in SB1 and MB1 settings can be easily understood at least qualitatively. As we 
learned in Sec. VC by lacking ~ 1/E terms, reach to high-energy asymptotic behavior is 
relatively fast. At low energies, sin 2 26 I * I 3 I is constrained to be small as one can show by 
the formulas based on the matter perturbation theory given in Appendix |Xj The first or- 
der correction term, from which the energy dependence comes in is small, of the order of 
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FIG. 11: The energy dependence of sin 2 29™ (upper three panels) and sin 2 2#[g (lower three 
panels) are plotted for the three typical settings SB1 (left panel), MB1 (middle panel), and MB2 
(right panel) defined in Sec. VIA The true value of #13 is taken as sin 2 2^13 = 0.05, which is 



indicated by the horizontal solid line in the figure. 



A/A 31 ~ 0.06 - 0.07(0.2) for SB1 and MB2 (MB1) settings. Whereas for sin 2 29 n there is 
no small parameter which forces it small. The mild energy dependence and the pinning to a 
small value makes sin 2 29 m small in the entire region of E. The relatively fast reach to the 
asymptotic behavior can be seen in most of the plots of energy dependence of the degeneracy 
solutions for SB1 and MB1 settings. Notice, however, that the asymptotic energy can be 
reached at much higher energies for MB2 setting. 

Here are comments on the solution (stv,5iv) : Riv essentially looks like Ru apart from 
the presence of no-solution regions. Given smallness of Rm, Riy must looks like Ru because 
they are the intrinsic degeneracy pairs, as discussed in Sec. |V} It is also true that the 
energy dependence of sin 2 29™ is very similar to the behavior of sin 2 29y 3 . Therefore, lifting 
degeneracy between 9™ and 9™ can be done with spectrum analysis via a similar manner 
as in the case of intrinsic degeneracy. If powerful enough the spectrum informations would 
solve both the degeneracy between the true solution and 9^, and the one between 9™ and 



9™ at the same time. 



Next, we discuss and D^j which are presented in the upper and lower three panels, 



respectively, in Fig. [12} We note that they are small in SB1 setting, leaving the sign- 
Am^ degeneracy intact in this short baseline setting. Notice, however, that it is not all 

(2) 

bad, because the smallness of D u { implies that no severe confusion takes place between 
CP violation and CP conservation. Now, the difference between SB1 and MB1 settings 
further develops in particular in large #13 region. The clear distinction between SB1 and 
MB1 settings is also prominent in the energy dependence presented in Fig. 13 Of course, it 
is basically due to larger matter effect in MB1 setting. It is interesting to observe that the 
difference shows up first in 5, but not quite for 6*13 at large fli 3 . 

In the SB1 and MB1 settings, as can be seen in Fig. I12I and D^j are largest in 
region of the largest possible #13 for which the sign-Am 2 ^ degeneracy solution exist. In this 
region 5 ~ n/2. On the other hand, Rm and i?rv are small in the region as is seen in 
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III 



5 tmc )]/vT (upper three panels) and Z?jy 



^IV_ (5 truey 7r ( lower three 



panels) denned in (75) is presented in sin 26*i3 — 5/vr space for the three typical settings SB1 (left 



panel), MB1 (middle panel), and MB2 (right panel) defined in Sec. VIA The regions of white 
color denote the regions of no sign-degeneracy solution. 
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FIG. 13: The energy dependence of the ratios DtJ defined in (75) are plotted for the three typical 
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D§> is 



settings SB1 (left panel), MB1 (middle panel), and MB2 (right panel) defined in Sec. 
not shown because the similarity in its energy dependence, apart from reversing the positive and 
negative regions of the ordinate. The true value of #13 is taken as sin 2 2#i3 = 0.05. 



Fig. 10 It is easy to understand these features. At around the largest value of #13 which 



allows the sign-Am^ degeneracy the two ellipses, the ones with normal and inverted mass 



hierarchies, barely overlap with each other. See Fig. 8 in [22] • The general feature of the 
bi-probability plot |24J tells us that in the overlap regions of the two ellipses the point of 
5 ~ 71 / 2 in the positive Am^ ellipse is close to point of 5 ~ 3n/2 of the negative Am^ 
ellipse. Therefore, ~ 1 and £)jy ~ 1 hold, explaining the above features. Because the 
center of the two ellipses are located at almost the same distances from the origin (which is 
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equal to S2 3 sm 26*13), sin — srv ~ s\. 

We want to note that the energy dependences of D^j and Djy are quite mild in energy 
region above the first oscillation maximum for SB1 and MB1 settings. Because the energy 



dependence of Dry is similar to that of D\tJ in for SB1 and MB1 settings, apart from 



)(2) 



reversing the positive and negative regions of ordinate, we do not present it in Fig. 13 
(Hereafter, we just quote the reference either [82] or [79] if the omitted figures are available 
in them.) For MB2 setting qualitative features of rapid up and down are very similar in 
all the figures of energy dependences presented, or abbreviated^ Considering the almost no 

and noting that spectrum 



energy dependence of sin 2 29 m 



sin" 



as given in Fig. 
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analysis is highly challenging at low energies, it would be difficult to resolve the sign-Am^ 
degeneracy by a single detector setting of either SB1 or MB1. 

We notice that the difference between SB 1-MB 1 and MB2 settings is always evident as 
can be seen in Figs. 



10, 11, 12 and 13 



Therefore, MB2 setting alone may have chance to 
resolve the sign-Am 2 ^ degeneracy [531 I57)] . Or, if the informations gained at around the 
second oscillation maximum can somehow be combined it would greatly help resolving the 
sign- Amjj! degeneracy [52TI541 157] . It may be expected even from our formulas obtained for 
a "mono-energetic neutrino beam" because the parameter regions with degeneracy solutions 
in SB1 and MB2 settings tend to "repel" (avoid to overlap) with each other at large #13, 
though not completely. 

Finally, we should note that abrupt termination of lines in the figures that appears in 
Figs. [IT] and [13] are either due to disappearance of the degeneracy solutions, or switching 



phenomenon between solutions that takes place due to our convention of labeling degeneracy 



solutions. The feature will be seen also in the foregoing subsections. See Sees. Ill D and IV 
for discussion on this point. 



D. Intrinsic and sign-Am 2 ^ degeneracies in the false 6*23 octant 

Now, we turn to the 823 octant degeneracy with solutions which lives in the different 6*23 
octant from the true solution. Having the overview at hand, we present the intrinsic and 



the sign-Am^ degeneracy solutions at the same time. Presented in Fig. 14 in sin 2 29\z — 5/ti 



space are the ratios Ry = (sin 2 29^ 3 — sin 2 291™°) / sin 2 26^g ue (top three panels), Ryi (next to 



top panels), Ryu (next next to top panels), and -Ryin (bottom three panels) defined in (74) 



for three typical cases of energies and baselines, SB1 (left panels), MB1 (middle panels), and 



MB2 (right panels) defined in Sec. VIA In Fig. [15| the energy dependence of degeneracy 
solutions of sin 2 26% is presented. We note that the energy dependence of sin 2 29™ is quite 
similar to that of sin 26%. The relat ion between the energy dependences of sin 2 26% and 
sin 2 26% (sin 2 29™ and sin 2 26% m ) is similar to the one between sin 2 29™ and sin 2 29™ given 
in Fig. [TTJ (See [79l E2].) It is natural because they are the intrinsic degeneracy partners, 
and hence they are not shown. 



One of the most notable features in Fig. 14 is again quite distinct behaviors in the MB2 



setting. In general, i?N are large (apart from the strips where i?N switches its sign) with 
notable exceptions of Ry for SB1 and MB1 settings, and i?vn f° r SB1 setting. It is also 
notable that behavior of Ryi an d -Rvm is reminiscent of the one of Ru in Fig. [6] in SB1 and 
MB1 settings in small 6 ) 13 region, except for the presence of no-solution region. Considering 
the small values of Ry and i?vn (except for i? V n f° r MB1), it is quite natural to see the 
behavior given the fact that they are the intrinsic degeneracy partners. The feature that 
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Ryi and i?vm trace the behavior of their intrinsic degeneracy partners also applies to MB2 
setting. The behavior of Ry and Ryu (Ryi an d -Rvm) is somewhat similar, apart from the 
presence of no-solution region, to that of Ru for SB1 (MB2) setting. It may be understood 
by the similar consideration using the bi-probability plot. Therefore, we concentrate below 
on SB1 and MB1 settings. 

As mentioned above Ry for SB1 and MB1 settings, and i?vn for SB1 setting are small 
in region sin 2 2^13 > 1CT 2 . The region of #13, however, is nothing but a good target for 
superbeam experiments, and it will be a challenge for them to lift the degeneracy solutions. 
Here, we try to understand this feature, but in a wider perspective which includes the 
energy dependence of sin 2 26^ and sin 2 29Yi l f° r SB1 and MB1 settings. In Fig. 15, we 
observe that the difference between sin 2 26 , } 3 (or, sin 2 26 , } 3 11 ) and sin 2 26^g ue is nonzero but 
energy independent in a wide region except for at very low energies, a somewhat unexpected 
behavior to see. 

Now, we point out that the behaviors mentioned above can be understood by formulating 
the 623 perturbation theory, as done in Appendix [Bj Namely, one can derive the perturbative 
expression of sin 2 2#f 3 and other quantities by assuming that deviation of #23 from 7r/4 is 
small, 823 — 7r/4 = e oct <C 1. One can expect that the expansion by e oct is indeed a good 
approximation because e.g., e oct = 0.05 for 623 = 42°. Then, we obtain 

sin 2 29^3 - sin 2 20f 3 ue = 4e oct sin 2 2#f 3 ue (76) 



for which we have used the fact that the last term in (B2) is negligibly small as far as we 
remain in a region sin 2 26^ c 3> Z. In fact, we confirmed that the correction term becomes 
non-negligible at small #13 around sin20*3 Ue = 10~ 3 . Certainly, the condition is fulfilled for 
the settings SB1 and MB1. The similar equation holds for but with replacement of 6\™ Q 
to ^{g 1 because their relation as the intrinsic degeneracy partner. Therefore, Ry and Ryu are 



small and the difference in (76) is approximately energy independent. One may ask why the 
feature does not exist in MB2 setting with small #13. As mentioned before, the solar term 
is dominant in this region. Therefore, if e oc t^± is comparable to s± the difference sy — si is 
no more small. 

Another notable point is that region of no degeneracy solution is not additive, as can 



be seen by comparing Figs. [10] and [14} That is, the region of no degeneracy solution with 
Amf^-sign and 623 octant flips (VII and VIII) is not the union of no-solution regions of the 
sign-Am^ (Ill) and the 623 octant (V) degeneracies. It is simply because the degeneracy 
solution with both sign and octant flips can exist even in a region of # 13 and S where e.g., 
the octant degeneracy solution does not exist. 



In Fig. 16, from the top to the bottom, the normalized differences between the true and 
fake solutions of phases, Dy \ Dyl, -Dyii; an d -^viii defined in (75) are presented for SB1 



(left panels), MB1 (middle panels), MB2 (right panels) settings. In Fig. 17, the energy 
dependence of Dy^ is plotted. The behavior of Dyl is similar to given in Fig. 9 , while 
those of Dy\ and Dy\ n are very similar to the ones of D^l and D^y (the latter not shown 
but the behavior explained) in Fig. 13 See [791 182]. We notice that for SB1 setting D$ (i 



either 1 or 2) is small in most of the regions of true values of 5 for all the solutions V-VIII. 
For MB1 setting the same statement applies for the solutions V and VI. A notable feature 
is that Dy^ (and -Dyn) * s sma U i n MB2 setting. It can also be understood from the 623 
perturbative formula for 6 given in Appendix [Bj The difference between 8\ and Sy is always 
suppressed by e oct . 
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FIG. 14: Plotted are the ratios Ry = [sin 2 2(9^ - sirr 



20f 3 ue 



]/ sin' 



26\f e 



(top three panels), Ryi 



(next to top panels 
frl in sin 2 26 



'13 



i?vn (next next to top panels), and -Rvm (bottom three panels) defined in 
5/tt space is presented for three typical cases of energies and baselines, SB1 (left 
and MB2 



panels), MB1 (middle panels), and MB2 (right panels 
color denote the regions of no degeneracy solution. 



defined in Sec. VIA The regions of white 
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FIG. 15: The energy dependence of sin 2 20^ are plotted for the three typical settings SB1 (left 



panel), MB1 (middle panel), and MB2 (right panel) defined in Sec. VIA For the behavior of the 
other solutions not shown, sin 20$, sin 2 20™ an d sin 2 20™ 

see the text. The true value of #13 
is taken as sin 2 20i3 = 0.05, which is indicated by the horizontal solid line in the figure. The true 
value of (9 23 is 42°. 



As discussed above the energy dependence is very mild for most of the solutions V— VIII, 
except for at low energies, E % 0.4 GeV, in SB1 and MB1 settings. Therefore, it may be 
extremely challenging for experiments with the settings to lift the degeneracy. Because of the 
likely difficulty in resolving the 2 3 octant degeneracy several methods have been proposed; 
the reactor-accelerator combined method [HI EH], the various ways to detect solar Am^ scale 
oscillations, using atmospheric [SnHH2] or accelerator neutrinos [MHSE], or both combined 
[oTj lrJ3"] . The silver channel could be of help |64j. As in the previous cases the behavior of 
degeneracy solutions are far more violent in MB2 setting. It by itself might mean the great 
sensitivity to resolve the degeneracy. Or, it is a natural way of thinking to combine it with 
the measurement at the first oscillation maximum. 

We give here a brief summary of the characteristic features of the degeneracy in super- 
beams, SB1, MB1, and MB2 settings. 

• A prominent difference between the true and the clone solutions at relatively large 
#13, sin 2 20 13 > 10 -2 , exists in in the intrinsic degeneracy, while for the sign-Am^ 
degeneracy it is in the phase difference D$ for SB1 and MB1 settings. 

• The solutions III for the sign-Am 2 ^ degeneracy appears to be difficult to resolve for 
SB1 setting even if spectrum information is available, because energy dependences are 
so weak for both Rm and -D m • The similar difficulty exists for MB1 if # 13 is large, 
sin 2 2#i3 > 1CF 2 . For the same reason, the solutions V and VII of the #23 octant 
degeneracy is difficult to lift. 

• The short baseline SB1 option is unique among the three superbeam settings in the 
sense that it by itself may not be able to lift the sign-Am^ and the #23 octant degen- 
eracies, but can provide a clean discovery of CP violation without confusion with CP 
conservation. This is in accord with the basic motivation for low energy superbeam 

[BB]. 

• In comparison with SB1 and MB1 settings, the features of degeneracy solutions are 
always quite distinct at MB2 setting, where the energy region around the second 
oscillation maximum is explored. It by itself, or combined with other settings, would 
provide ways to help resolving the eightfold degeneracy. 
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FIG. 16: From top panels to bottom panels presented are in order: Dy = (5 V — <5 truc )/7r, 
L>yj = [5 V1 - (tt- 5 tmc )]/vr, D { yl = [5 Yn - (vr - <5 true )]/vr, and D$ n = (5 vm - <5 true )/^ denned in 
(75) for three typical cases of energies and baselines in sin 2 26\™ e — 5 trnc /7r space. The regions of 
white color denote the regions of no sign-degeneracy solution. 
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FIG. 17: The energy dependence of the ratio 1)1 defined in (|75l) are plotted for the three typical 



settings SB1 (left panel), MB1 (middle panel), and MB2 (right panel) defined in Sec. VIA For 

shavior of the other 
of is taken as sin 2 2#i3 = 0.05. 



the behavior of the other solutions not shown, D yj , and D yjjj, see the text. The true value 



E. Parameter degeneracy in neutrino factory setting 

In this subsection, we display the features of various degeneracy solutions by taking 
a setting which may be appropriate for neutrino factory. The two-detector setting with 
baselines L = 3000 — 4000 km and L ~ 7000 km seems to be considered as the "standard" 
one |83J both for measurement of standard mixing parameters [231 EI], possibly as well as 
for search for effects of NSI [3H ESI EE]- For the former purpose, the far detector at the 
"magic baseline" (as named by [HI]) plays a key role in resolving the degeneracy because of 
absence of 5 dependence [27]. By restricting our purpose to illuminate the features of the 
degeneracy, we use the setup with just one detector at L = 4000 km. 

Note that we have used set of probabilities Pivp, — > v e ) and Pip^ — > v e ) to obtain the 
degenerate solutions. Therefore, if you want to consider the more realistic setting of a 
neutrino factory in which T-conjugate (golden) channels will be used, please regard 5 as 
2ti-5. 



In Figs. 18 and 19, the differences between the true solution and the fake degeneracy 
solutions, i?N and (N=II— VIII, i = 1 or 2), respectively, are plotted. As can be seen 
in these figures the differences between the true solution and the fake degeneracy solutions 
are generically larger than the cases of superbeam type settings discussed in the previous 
subsections. In accord with the expected higher sensitivities, we extend the region of sin 2 2#i3 
to 10" 4 . 
In Figs. 



sented for t 



20 



and 



21 



the energy dependences of sin^ 20f 3 and D$, respectively, are pre 



; iree representative solutions, II, III, and V. For energy dependences of the other 
solutions see [7JEE2]- Generally speaking, the energy dependences of both of the quantities 
are significant compared to those in the SB1 and MB1 settings. Notable exceptions are the 
solution V (both sin 2 29Y 3 and Dy ), and possibly sin 2 2#J I 3 I and sin 2 26^ n , all except for the 
low energy region E ^ 10 GeV. Unless there is a sensitivity to the low-energy region it 
would be difficult to resolve the degeneracy, in particular V, by the spectrum informations. 
Therefore, it is extremely important to lower the threshold into E < 10 GeV to resolve the 
degeneracy by spectrum analysis. An extensive effort toward this direction is made and the 
task is in progress [57] . 

It is possible to understand, at least qualitatively, smallness of the difference from the 
true solution and lack of strong energy dependence of the solutions V and VII. Because the 
623 perturbation theory also applies to NF setting, it can be expected that the difference 
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FIG. 18: The ratio i?n— -Rvill are plotted in sin 2 26*13— <V 7r space for a typical baseline and energy of 
neutrino factory setting. The ratio i?N is defined as i?N = [sin 2 28^ — sin 2 29\™ c ]/ sin 2 29\™ e defined 
in d74b 
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FIG. 19: The normalized differences , -Djfj, -D{v> > ^Vl > ^vn> ana - ^vin are plotted in 
order in sin 2 2#i3 — 5/ir space for a typical baseline and energy of neutrino factory setting. Z?n is 
defined in J75I as = [5 N 



(vr - 5 tmc )]/vr and D$ = [8 N - 8 tTUC }/7r. 



between the true solution and the clone one V is small. The similar statement holds for 
the solution VII given the smallness of the energy dependent term in sin 20} I 3 I . Then, the 
question is why sin 2 26*™ is small and lacks the significant energy dependence despite that 
the matter perturbation theory is not valid for NF setting. Qualitatively, the answer is that 
pinning to a small value due to the fact that the assumed true value itself is small, and lack 



of energy dependence because of fast reach to the asymptotic behavior discussed in Sec. VC 
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FIG. 20: The energy dependences of sin 2 20^ (left panel), sin 2 (middle panel), and sin 2 20^ 
(right panel) are plotted. The true value of #13 is taken as sin 2 2#i3 = 0.001, which is indicated by 
the horizontal solid line in the figure. 
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FIG. 21: The energy dependences of (left panel), (middle panel), and Dy^ (right panel) 



are plotted. D$ and D-jp are defined in (|75|). 



One notices that the intrinsic solution II in NF setting has the similar features as the one 
in MB2 setting, as can be seen by comparing Figs. 18 and 19 to Figs. |5]and[6j It is because 
the value of #13 taken is small and the atmospheric term is comparable to the solar term. 
Parallelism is not so complete in the other types of degeneracies, but some features can be 
understood in analogy to the case of MB2 setting. 14 In doing so the difference due to the 
much wider no-solution region of the sign-Am^ degeneracy due to long baseline must, of 
course, be taken into account. 



VII. PARAMETER DEGENERACY IN T-CONJUGATE MEASUREMENT 



We analyze in this section the problem of parameter degeneracy in T-conjugate mea- 
surement in neutrino oscillation. Though measurement of T violation does not appear to be 
feasible immediately understanding its structure may be interesting theoretically. The topics 
was first treated in [29], but we make the structure of the degeneracy more transparent in 
harmony with the symmetry argument in Sec. II B In fact, the structure of the degeneracy 
with T-conjugate measurement is one of the simplest one among the cases discussed due to 



An example is that there is the region that Ry is large in small #13. This is for the same reason of 
the case in MB2 setting, the difference between the true and V of order ~ e ct^±j is negligible compare 
with si. Another example is much stronger energy dependence of su in NF setting than SBl's, which is 
reminiscent of the feature of MB2 setting. It comes from larger effect of the solar-scale oscillation term. 
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the symmetry. It is also an ideal tool to obtain the vacuum limit. 



A. The intrinsic degeneracy in T-conjugate measurement 



The intrinsic degeneracy solutions (sj, Si) (i=l, 2) are defined in —> v e channel by 
P — Z = X±s\ + Y±s\ (cos 8\ cos A 3i =)= sin 8\ sin A 3 i) , 



P — Z = X±s 2 + Y±s 2 (cos <5 2 cos A 31 =f sin 8 2 sin A 31 ) 

and in T-conjugate channel by 

P T — Z = X±s\ + Y±s\ (cos^i cos A 3 i ± sin^i sin A 3 i) 
' = X±sl + Y±S2 (cos5 2 cos A31 ± sin<5 2 sin A 3i ) 



P 



T 



(77) 



(78) 



By subtracting two equations in (77) and (78), respectively, and then subtracting and adding 
the resultant two equations, we obtain (assuming Y± sin A 31 ^ 0) 

s 2 sin<5 2 = Sisin^i, 



s 2 cos 82 = Si cos 5i + 



X 4 



Y+ cos A 



(79) 
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Inserting (79) into cos <5 2 + sin 5 2 = 1 yields the equation for s 2 in a form (s 2 — Si)(s 2 — s 2 j) 
0, which admits the intrinsic degeneracy solution 



sj + 2 



Y± cos A 



31 



X. 



± 



si cos 8\ + 



Y± cos A 



31 



x 4 



(80) 



Given the solution s 2 = sn the solution for <5 2 can be obtained by using ( 79 ) as 

snsin^n = Sisin^i, 



sn cos 8n — — ( si cos 81 + 



Y± cos A 3 i 
X + 



(81) 



By further expanding by assuming it small, the solution obtained in [2S] is reproduced; 



l±cosA 31 
sn — «i H v cos 8 X . 

A + 



(82) 



B. The sign-Am 2 degeneracy in T-conjugate measurement 



As we learned in Sec. II B the symmetry argument tells us that there exists the sign- 
Am 2 degeneracy in T-conjugate measurement. In this subsection we verify it by deriving 
explicit solutions without recourse to the symmetry argument. We denote Si 3 variable for 
the opposite-sign Am 2 x solution as (s 3 , 8 3 ), whose two (as we prove) solutions will be denoted 
as (sin, (5in) and (srv,<$rv)- 

The sign-Am^ degeneracy is defined by the following two sets of equations: 

P — Z = X±s\ + Y±si (cos 81 cos A 31 =1= sin 81 sin A 31 ) , 

P — Z = X T sl + KpS 3 (cos<5 3 cos A 31 ± sin 5 3 sin A 31 ) , (83) 
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and 



P T — Z = X±s\ + Y±si (cos 5i cos A 3 i ± sin 5i sin A 31 ) 
P T — Z = X T s\ + Y T S3 (cos S3 cos A31 =F sin £3 sin A31) 



(84) 



By the similar procedure as in the previous subsection we obtain (assuming Y± sin A 31 7^ 0) 



s 3 sin S 3 
s 3 cos 5 3 



Si sin 5 



i- 



si cos<5i 



1 



Y T cos A 3 i 



(X ± sl - X T s 2 3 ) 



(85) 



Inserting (85) into cos S3 + sin <5 3 = 1 leads to the equation for S3 as 



2 _ X± 2 



2 _ x± 2 

53 x T Su 



0. 



(86) 



where sn is defined in (80). We note that the relation ^ is essential to reduce the equation 
in (86) to the current form. The solutions for CP phase 5 can be obtained by inserting 
the solutions of (86), its positive root, into (85) with use of ([9]), and for S\y together with 
(79). The solutions of 5 together with the 6*13 solutions of the sign-Am^ degeneracy are 
summarized as below: 



srv 



IX* 

x^ 



IX 



± 



S m = n-5 1 , 

-.s n . 8jy = TV — 5ji- 



(87) 



The solutions (87) are in perfect agreement with the expectation based on invariance of the 



oscillation probability given in Sec. II B Figure [2] clearly exhibits the structure obtained in 
(|87|). 



C. The octant degeneracy in T-conjugate measurement 

The 6*23 octant degeneracy in T-conjugate measurement is defined by 



P 
P 



x^true 2 



7truc 



P 



*7 t- i±si (cos 5i cos A 3 i =F sin Si sin A 31 ) + Z 

Y±s 5 (cos S 5 cos A 3 i =F sin S 5 sin A 3i ) + Z false , 
X^ uc s'i + Y± Sl (cos ft cos A 3 i ± sin<5i sin A 31 ) + Z true , 



false J2 
± S 5 
T vtrue 2 



p i = x^s* + Y±s 5 (cos S 5 cos A 3 i ± sin S 5 sin A 3i ) 



7 false 



From these equations, we obtain 

S5COS55 = sicos5i + 
s 5 sin5 5 = Sisin^i. 



1 



Y± cos A 3 i 



true „2 
'1 



x^false 2 i rytrue 
A ± S 5 + Zj 



7 false! 



(88) 
(89) 



(90) 
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The relation cos 2 85 + sin 2 <5 5 = 1 gives the quadratic equation for s\ which results in the 
similar solution 



'v,vi 



(x4 alse ) 2 



[Xf se {X^sl + Y± cos A31S1 cos St + {Z l 



^falsc 



)} 



+ -(y ± cosA 31 ) 2 t 4 



^T-oct-intr 
oct-intr\2 



(91) 



where the upper — (lower +) sign is for sy (svi)- The functions £)T-oct-mtr anc i ^T-oct-mtr are 
defined by 



^T-oct-intr = (y ±cos A 31 ) 2 



X 



(x^ he Sl cos Si + cos A 3 i^ + X£ alsc I (Z tTUC - Z falsc ) - (X^ uc - x£ alsc )s 2 J 



^T-oct-intr = ^ ^T-oct-intr 
e 2 3-S>7r/4 



Y± cos A 31 ( X±Si cos 5i + -Y± cos A 31 



(92) 



where X± and Y± in the last line is meant to be those at 623 = it /A. Once the solutions sy 
and syi are known one can readily obtain Sy and Syi by inserting the S13 solutions into (90). 



Though, we do not discuss in any detail, it must be obvious that the solutions (syn, <5vn) 
and (syiii, 6Vm)> the ones with both the Am^-sign and the octant flips, are given by using 
the solutions obtained in this subsection by the same type of equation as (87): Syn = 

Sy, 



^J^falsc ^J^falsc 



Sy, Syiii 



7T - 



and Syin = TT — Syi. Alternatively, 



they can be also obtained by the octant flip mapping, implicitly given in (91) and (90), from 
(■sin, 5m) and (srv, ^rv) obtained in the previous section. 



D. Parameter Degeneracy in Vacuum 

Though it is pedagogically useful to work out the parameter degeneracy in vacuum we 
just give the results by taking the vacuum limit in the degeneracy solutions obtained for T 
conjugate measurement. We need the results in vacuum to define our convention we used 
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to specify unambiguously the solutions for the sign-Am^ degeneracy in Sec. Ill D 

In fact, it is straightforward to observe that the expressions of the degeneracy solutions 
in vacuum: They are almost identical to those obtained in this section; The only necessary 
step is to take the vacuum limit 

limX ± = X vac = 4s 2 3 sin 2 A 3 i, 

a— >0 

lim Y± = ±Y vac = ±2 sin 26\2 sin 2#2 3 A2i sin A 3 i, 

a— >0 

limZ = Z vac = sin 2 26*12, (93) 

a— >0 



15 Notice that the mass hierarchy does matter in the discussion of parameter degeneracy in vacuum though 
it might be thought contrary. In fact, it is known that the hierarchy can be determined by measuring the 
sign of the solar-atmospheric interference term in vacuum [88] . 
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in the solutions (80), (81), and (87). Their explicit forms are given in Appendix D of [38J. 



The degeneracy solutions which involve #23 octant flip can also be obtained by taking the 
same limit in the solutions obtained in Sec. IVII CI 



VIII. PARAMETER DEGENERACY WITH GOLDEN AND SILVER CHANNELS 

We discuss the parameter degeneracy for a given measurement in the u e — > (golden) 
and the v e — > v T (silver) channels. The oscillation probability in the former and the latter 
channels are given by P T in (pi) and P s in (10), respectively. We note that detection of 



v T requires the energy at least above r production threshold, and therefore most probably, 
neutrino factory would be the appropriate place for the silver channel [71] . 



A. The intrinsic degeneracy in Golden-Silver measurement 

The intrinsic degeneracy is defined by 

P T — Z = X±s\ + Y±si (cos 5i cos A 31 ± sin S 1 sin A 31 ) 
P T — Z = X±s\ + Y±S2 (cos £2 cos A31 ± sin5 2 sin A 3i ) 



(94) 



and 



P s - tan 2 6^Z 



723ZJ = cot 2 023X±sl — Y±si (cos Si cos A31 ± sin Si sin A 3i ) 
P s — tan 2 623Z = cot 2 #23Ar± s 2 — Y±s 2 (cos 5 2 cos A 31 ± sin 5 2 sin A 31 ) 



By subtracting two equations in (|94|) and (|95|) we obtain 
X ± (sl-sl)+Y ± Sj 

where 



0. 



X±(4~4) -tan 2 9 23 Y ± S I = 0, 



Sf± = [cos A 31 (si cos Si — S2 cos S 2 ) ± sin A 31 (si sin Si — s 2 sin <5 2 )] 



(95) 



(96) 



(97) 



Then, it follows that sn = si (using the label sn for the intrinsic solution) and Si — 
assuming that Y± 7^ 0. The former result is, of course, expected by the feature of "shrunk 



ellipse" in the bi-probability plot given in Fig. 22 Using su = Si the equation SV = has a 
solution, apart from the trivial solution 8% = 8%, as 

cos #11 = cos(<5i =F 2A 3 i), sin<5 n = — sin(5i =F 2A 31 ) (98) 

which implies that 

Sn = 2tt - (Si =F 2A31) (mod. 2tt). (99) 
The structure of the solutions of S should be obvious from the form of the oscillation prob- 



abilities in (94) and (95) which is reflected to the feature of shrunk ellipse in Fig. 22 



The two degenerate solutions must have the same values of cos(<5 =F A 31 ), and hence 
5nTA 31 = 27r-(5 lT A 31 ). 
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L=4000km, E=20GeV 
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FIG. 22: P T — P s (golden-silver) bi-probability plot for a neutrino factory setting. Black solid 
line : sin 2 2#i3 = 0.004(0.001) with normal hierarchy for larger (smaller) oscillation probability. 
Blue dashed line : sin 2 2#i3 = 0.0158(0.0039) with inverted hierarchy for lager (smaller) oscillation 
probability. Notice that black and blue lines are overlapping. #23 is taken to be 42 degree. 



B. The sign-Am 2 degeneracy in Golden-Silver measurement 

The sign-Am 2 degeneracy is denned by 



P T - Z 
P T - Z 



X±s\ + Y±si (cos Si cos A 31 ± sin Si sin A 31 ) 
X T sl + Y T S3 (cos #3 cos A 3 i =f sin 5 3 sin A 3i ) 



(100) 



P s — tan 2 6*23^ = cot 2 6* 2 3^± s i — Y±si (cos^i cosA 31 ± sin#i sinA 31 ) 
P s — tan 2 #23-^ = cot 2 ^23^=f s 3 ~~ ^t s 3 (cos £3 cos A31 =F sin £3 sin A31) 



By subtracting two equations in (100) and (101) we obtain 



X ± sl-X T sl + S & 



0. 



X±s\ - X T s\ - tan 2 9 23 S S = 0, 



where 



Then, (102) implies that the two separate factors both have to vanish: 



sin 



x^ 



-si, 



s. 



(101) 



(102) 



Ss = cos A31 (Y±si cos^i — F T S3 cos £3) ± sin A31 (Kj-Si sin^i + KpSssin^) . (103) 



(104) 



where we use the notations (sm, Sm) and (srv,<5rv) for the sign-Am^ degeneracy solutions. 



The former result is in agreement with the symmetry argument given in Sec. |II B[ Because 
sii = si, siv = sm holds. 

By using the above S3 solution and the relation (|9| in Se c. \H\ it is easy to show that the 
equation Ss = can be converted to the same equation as (98) (after obvious replacement 
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of 82 to 83) apart from the sign change of cos £3 term. It means that 7r 
the same equation as (99). Then, we obtain the solutions 



S3 obeys exactly 



£111 



7T 
7T 



Si, 

811 



-it + (81 =F 2A31) (mod. 2?r). (105) 
The structure of the solution in one-to-one correspondence to the intrinsic degeneracy solu- 



tions is perfectly consistent with the symmetry argument in Sec. II B 



C. The octant degeneracy in Golden-Silver measurement 



Now, we discuss the #23 octant degeneracy. We will see that simplicity of the golden-silver 
setting prevails in it. The #23 octant degeneracy solutions (s 5 , 8 5 ) satisfy 



P T = X^ ue si + Y±S! (cos Si cos A 3 i ± sin S x sin A 31 ) + Z 



true 



P 1 = X± lse s 2 5 + Y±s 5 (cos 8 5 cos A 3 i ± sin 8 5 sin A 3 i) + Z 



false 



(106) 



and 



p s 



cot 2 ^23^± ue -5i - Y±si (cos 81 cos A 31 ± sin 81 sin A 31 ) 
tan 2 2 3-X± lse Sg - Y±s 5 (cos 8 5 cos A 31 ± sin 5 5 sin A 3i ) 



tan 2 e 23 Z truc , 
- cot 2 6 23 Z Msc . 



(107) 



Using (39), (107) can be written as 

p s = X i^e s 2 _ Y±8l (cqs^ cos A 31 ± sin ^ sin A 31 ) + Z false , 



P" 



X± ne sl - Y±s 5 (cos 8 5 cos A 3 i ± sin 8 2 sin A 3i ) 



?true 



(108) 



By subtracting two equations in (106) and (108) with the same octant labels one can easily 
obtain the solution of S5 as s 2 = sf which leads to 



sv = Syi = Si. 

To obtain 8$ , we proceed as usual which leads to the result 

1 



(109) 



cos(cV T A 3 i) = 008(5! =F A 3 i) + — [{X^ - Xi alse ) s 2 + (Z true - Z false )] . (110) 



Sl Y 

As dictated by the general argument, syi must be given by sy as above, and Syi as 

Syj = 2tt - (8 V =F 2 A31) (mod. 2tt), 



111) 



using (99) because they are the intrinsic degeneracy pair. In fact, it is easy to see that Sy 
and Syi are the two solutions which satisfy cos(cV =F A 31 ) = cos(cVi T A 31 ) in (110). 

As in the case of T-conjugate measurement described in Sec. VII C the solutions 
(svii,5vii) and (svn, Syn), the ones with octant as well as the Amf^-sign flips, can be ob- 
tained from the above ones as 



Syn 
Syn 



s VIII 



/ J^false 
false 



sv, 



7T 



'VIII 



7T 



?VI- 



(112) 
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IX. PARAMETER DEGENERACY IN CPT-CONJUGATE MEASUREMENT 



We discuss in this section the problem of parameter degeneracy in CPT violation measure- 
ment. CPT- violation observable was considered to be useful to resolve the mass hierarchy 
because it gives the probability difference P — P CPT which is largest among the similar 
quantities [89| 190]. 



A. The intrinsic degeneracy in CPT-conjugate measurement 



With expression of the oscillation probabilities in (|3j) and (|6j), the intrinsic degeneracy 
solutions (sj, Si) (i=l, 2) in CPT-conjugate measurement are defined with — > v e channel 
by 



P 
P 



Z = X±s\ + Y±si (cos<5i cos A31 =1= sin<5i sin A31) 
Z = X±s\ + Y±s 2 (cos5 2 cos A 3 i =F sin S 2 sin A 31 ) 



(113) 



and in CPT-conjugate u e — > v„ channel by 



P PT — Z = X T s\ — Y T s\ (cos<5i cos A31 =F sin^i sin A31) 
P PT — Z = X T s\ — Y T S2 (cos S2 cos A31 =F sin5 2 sin A31) 



(114) 



By subtracting two equations in (113) and (114) respectively, and subtracting and adding 
the resultant two equations we obtain, assuming that 7^ 0, 



S 



it 



0. 



2{s\ - si) ± C ( - ) S_ 



it 



0. 



(115) 



where C^' and Sj± are defined in (17) and (97), respectively. We then obtain su = s 1; that 
is, the intrinsic degeneracy solution of #13 for CPT conjugate measurement is identical to 



the true one, in agreement with the expectation of the bi-probability plot. See Fig. [23] It is 

(116) 



obvious that the solution of the first equation (115) is given by 

(mod. 271"). 



5u = 2tt - (Si ± 2A 3 i) 



B. The sign-Am 2 degeneracy in CPT-conjugate measurement 

We now discuss the sign-Am^ degeneracy in CPT-conjugate measurement. The true 
input solution (si, Si) and the opposite Am^-sign clone solution (s 3 , 5 3) satisfy the following 
equations. In the — > v e channel, 

P — Z = X±s\ + Y±si (cos Si cos A 3 i =f sin Si sin A 31 ) , 

P — Z = X T sl + Y T s 3 (cos<5 3 cos A 31 ± sin 5 3 sin A 31 ) , (H7) 

and in CPT-conjugate channel 

pCPT _ g _ x T s\ — Y T s\ (cos Si cos A31 =1= sin Si sin A31) , 

p cpt _ z = x ± s2-F ± s3(cos53CosA 31 ±sin53sinA3i). (118) 
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L=300km, E=0.75GeV 
0.010, 




FIG. 23: P — P CPT bi-probability plot. The black and the blue solid (dashed) lines, which 
correspond respectively to the normal and the inverted hierarchies, are for sin 2 29\3 = 0.01(0.001). 
The region with small probabilities P < 0.0007 is magnified into the sub-panel to show more clearly 
the crossing of the two shrunk ellipses in the region. 



By subtracting two equations in (117) and (118), respectively, and then subtracting and 
adding the resultant two equations, we obtain 

— S3 cos £3 cos A31 =F S3 sin 63 sin A31 = 0, 



rpCPT 

2 i± 



CM 

T CPT _ Q„ 2 

— ± s 3 cos 5 3 cos A 3 i + s 3 sin5 3 sin A 31 = 0. 



c(-) 



(119) 
(120) 



where we have denned 



Ti± T = ± J E ( " ) s 2 + D (+) [sicos5icosA 3 i^sisin5 1 sinA 3 i] 



T. 



CPT 
2± 



E^ + 's\ ± ' [s\ cos 5i cos A31 =F si sin 5\ sin A 



31 



(121) 



where C^> is defined in (17), while and E^' are given in (27). From (119) and (120) 
it is straightforward to obtain the S3 solution: 



sin = siv 




[cm 



(122) 



where the ± sign is the hierarchy sign. Upon obtaining the S3 solution one can readily obtain 
S3 by solving (119) for cos(<5 3 =F A 31 ). The solutions read 



#111 
<5iv 



±A3i + arccos 



(. 



CPT 
1± 



VC (+) siii 



(mod. 2n), 



2tt - S m ± 2A 



*31- 



(123) 



One might have suspected, from the feature of the bi-probability plot in Fig. [23j that 
the sign-Am^ degeneracy solutions exist in a very limited region of small #13. Therefore, 
we present in Fig. 24 the region of no solution of the sign-Am^ degeneracy region by the 
shaded region. 
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FIG. 24: Depicted as the shaded areas in the sin 2 2^13 — 5/tt space are the regions where no 
sign-Am^ degeneracy solution exists for CPT conjugate measurement for MB1 (left panel) and 
MB2 (right panel) settings. The true mass hierarchy is taken to be the normal one. 



C. The 823 octant degeneracy in CPT-conjugate measurement 

The #23 octant degeneracy is denned by the following two sets of equations: 

P = X^ uc sl + KfcSx (cos 5 X cos A 3 i =F sin 8 l sin A 31 ) + Z tmc , 

P = X^sl + Y±s 5 (cos 8 5 cos A31 =F sin 8 5 sin A 31 ) + Z false . (124) 



pCPT = x true s 2 _ y^ ^ ^ CQg ^ ^ ^ ^ ^ ^ + ^truc^ 
pCPT = X false s 2 _ y^ ^ ^ ^ ^ ^ ^ ^ ^ ^ + ^falsc^ 



(125) 



Following the similar procedure as before it is not difficult to obtain the equation which 
involve neither 85 nor 81. Then, by using ^ and (39) we obtain 



Sy = svi = tan 6 



Then, the phase £5 is determined as 




^truo ^falsc 



/j^true J£"true 



(126) 



Sy C0S(5 V ± A31) = Sl CGS(5 1 ± A31) + - ^) ( ZtmC - ZfalSe ) • 



(127) 



As in the previous section the intrinsic degeneracy partner 8yi is given by using (116) as 

8 Y1 = 2tt - (8 V ± 2A 3 i) (mod. 2tt). (128) 

As in the previous cases, the solutions with octant as well as the Am^-sign flips, 
(svn, 8yu) and (svii,^vii), are given by the general argument as 



svn — S VIII — S>± 
Syn = 1]2 PT (sy,8y) 



£ PT (sv,6y) 



(129) 



where £± PT and ?7^ PT are defined in (122) and the first line in (123), respectively, as a 
function of (si, 81). 
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X. CONCLUSION 



In this paper, we have analyzed the problem of parameter degeneracy in various settings, 
CP-conjugate, T-conjugate, CPT-conjugate measurement, as well as combining the golden 
and the silver channels. Using the approximate form of the oscillation probabilities obtained 
by Cervera et al. we have derived, for the first time except for CP-conjugate setting, the 
exact analytic expressions of the eightfold degeneracy solutions in all these cases assuming 
#23 7^ 7r/4. We hope that the simple explicit expressions of the clone solutions nicely fill the 
"hole" of informations and help understand the nature of the degeneracy. Furthermore, they 
would prove to be useful if they can be implemented in an analysis codes such as [HI [92] 
to facilitate the search for fake minima of the x 2 . All in all, we expect that such solutions 
would help in correctly interpreting data to be taken in precision measurement in the future 
neutrino oscillation experiments. 

We have presented a new view of the parameter degeneracy as invariance under the 
discrete mappings of the flavor mixing parameters including the mass hierarchies. The 
explicit forms of the mappings can be obtained by the symmetry respected by a pair of the 
oscillation probabilities alone, (11) in Sec. II B, for the sign-Am^ degeneracy in T-conjugate 
and the Golden-Silver measurement. In all the other cases, the explicit forms of the mappings 
are given by the analytic expressions of the degeneracy solutions. Nature of the degeneracy 
as the intrinsic degeneracy duplicated by the sign of Am^ and #23 octant is now given its 
precise meaning as the mapping relations between each pair of degeneracy solutions given 
in (65). The structure emerged, the one-to-one correspondence between the true and the 
degeneracy solutions which we call the solution network, is illustrated pictorially in Fig. |4} 
We have also clarified the relationships between the degeneracy solutions for the given true 



mass hierarchies, normal or inverted, in Sec. VB 



The explicit analytic expressions of the eightfold degeneracy solutions are used to make 
plots of the difference between the true and the degeneracy solutions to give an overview of 
the degeneracy. The features of the degeneracy solutions are so profound, making the true 
overview of the degeneracy extremely difficult. However, we believe that we have illuminated 
some of the significant features by taking the three superbeam type settings as well as the 
one akin to neutrino factory. In particular, the mild energy dependence of the difference 
between the true and the clone solutions indicate the robustness of the degeneracy against 
spectrum analysis. These plots illuminate which degeneracy is likely to be difficult to lift, 
hence it could be useful to design future experiments in preparation of the degeneracy to be 
met in the measurement. Such precise understanding of the parameter degeneracy would 
be a definitive help if future precision measurement could be contaminated by new effects 
outside of the standard three-flavor mixing of neutrinos. 
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Appendix A: Matter Perturbation Theory of Sign-Am^ Degeneracy 

Here, we present the approximate formulas of the sign-Am^ degeneracy solutions within 
the framework of matter perturbation theory [77| ITS] which assumes <C l. 16 As can be 
seen in ([2]), the values of the ratio (assuming the matter density and Am^ referenced in 
the equation) are 0.060, 0.066, and 0.20, respectively, for SB1, MB2, MB1 settings discussed 
in Sec. |VI[ Therefore, the condition for validity of matter perturbation theory holds in a 
good approximation for the former two settings. For MB1 setting the approximation may 
be modest but we may utilize it for a qualitative discussions. 

H aving the analytic expressions of the sign-Am 31 degeneracy solutions at hand (see 



Sec. 



Ill C ) , it is straightforward to expand it in terms of the small parameter 



We 



only present the results. To first order in we obtain the following expressions. For 9 



13, 



sin 2 29™ = sin 2 2^3 



1 + 



A 4 sin 5i ( A 3i cos A31 - sin A 3i ) (2si A vac cos A 3i ± V vac cos Si 



A 



■31 



sin 2 A 3 i(2siX vac cos Si ± V vac cos A 31 ) 



(Al) 



and for S, 
cos Sm 



( A \ [^^(AaicosAsi-sinAax) 

cosdi+4 - — r^— 

I A ~ ' sm 2 A 3 i 



x 



\^3l/ L sm ^31 J 

g 2 A vac + Z vac cos(^i - A 3i ) cos(^i + A 3 i) ± siY vac cos Si cos A 3 i 

-l-cV nncA. X97 nr>cA„, 



iSl^vac COS Si + 2Z vac COS A 31 



sin Sm = sin Si + 4 



.4 



A 



31 



cos <5i ( A31 cos A 3 i — sin A 3 i ) 



sin 2 A 



31 



gj^vac + Z vax , CQS(^1 - A31) CQSyji + A31) ± SiY vax . cos Si cos A 3 i 

±SiY Yac cos Si + 2Z vac cos A 3 i 



(A2) 



where X var etc. are defined in Sec. VII D 



Appendix B: Perturbation Theory of #23 Octant Degeneracy 



In order to understand features of #23 octant degeneracy it is useful to have a perturbative 
framework assuming that deviation of #23 from the maximal is small, 

7 + 6 °ct ( 



9 



23 



(Bl) 



By expanding the octant degeneracy solution derived in Sec. IV we obtain #13 to first order 
in e oct , as 

sin 2 20*3(1 + 4e oct ) 



sin 2 2^ 



+16e oct Z 



±2g 2 v /XlA;sin2A3i + s x {Y ± sin(6i ± A 3 i) + Y T sin(^ =p A 3 i)} 
y/X±X^ [si{Y± sm(Si T A31) + Y T sm(Si ± A 31 )} =p 2Z sin 2A 31 ] 



(B2) 



16 It is known that the matter perturbation theory treatment of the parameter degeneracy gives rise to a 
transparent view of the degeneracy, which include e.g., decoupling between degeneracies 15H 155] . 
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where X± etc. implies those evaluated at 823 = vr/4. Similarly, we obtain for S 

cos Sy = cos 5\ + 2e oc t sin S\ 

x ~ {s\X T - Z)Y ± cos(5 1 ± A 3 i) + (4g± ~ cos(^i =F A31) - 2 Sl Z(X ± - X T ) ' 
s lv ^lX^{ Sl F T sin(5 1 ± A31) + SiKt sin(5 1 T A 31 ) - 2Z sin 2A 31 } J ' 

sin 5y = sin Si — 2e oct cos Si 

x " (^X T - Z)y ± cos(£i ± A31) + - Z)Y T cos(S 1 T A31) - 2siZ(X ± - X T ) ' 

X siy/X±X T {siY T sin(5i ± A 31 ) + siY± sin{S 1 =F A 31 ) - 2Zsin2A 31 } 

(B3) 
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